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ABSTRACT 

Aims. The X-ray-TeV connection and the evolution of the emitting particle population is studied in high-energy peaked BL Lac objects, by 
obtaining spectral information in both bands on sub-hour timescales. 

Methods. Simultaneous observations with HESS, Chandra and the Bronberg optical observatory were performed on the BL Lac object PKS 
2155-304 in the night of July 29-30 2006, when the source underwent a major 7-ray outburst during its high-activity state of Summer 2006. 
This event took place about 44 hours after the other major outburst of the night of July 27-28, which is known for its ultrafast variability. An 
unprecedented 6 to 8 hours of simultaneous, uninterrupted coverage was achieved, with spectra and light curves measured down to 7 and 2 minute 
timescales, respectively. 

Results. The source exhibited one major flare along the night, at high energies. The 7-ray flux reached a maximum of 1 1 times the Crab flux 
(>400 GeV), with rise/decay timescales of ~ 1 hour, plus a few smaller-amplitude flares superimposed on the decaying phase. The emission in 
the X-ray and VHE 7-ray bands is strongly correlated, with no evidence of lags. The spectra also evolve with similar patterns, and are always soft 
(photon index r>2), indicating no strong shift of the peaks in the spectral energy distribution towards higher energies. Only at the flare maximum 
is there evidence that the 7-ray peak is inside the observed passband, at ~400-600 GeV. The VHE spectrum shows a curvature that is variable 
with time and stronger at higher fluxes. The huge VHE variations (~22x) are only accompanied by small-amplitude X-ray and optical variations 
(factor 2 and 15% respectively). The source has shown for the first time in a high-energy peaked BL Lac object a large Compton dominance 
(Lc/Ls ~10) - rapidly evolving - and a cubic relation between VHE and X-ray flux variations, during a decaying phase. These results challenge 
the common scenarios for the TeV-blazar emission. 

Key words. Galaxies: active - BL Lacertae objects: Individual: PKS 2155-304 - Gamma rays: observations - X-rays: galaxies 



1. Introduction 



Send offprint requests to: luigi.costamante@stanford.edu or Among blazars, high-ene rgy peaked BL Lac objects (HBL; 
Rolf.Buehler@mpi-hd.mpg.de iGiommi & Padovanill 19941) are characterized by the highest en- 

* supported by CAPES Foundation, Ministry of Education of Brazil ergy particles and most violent acceleration processes. In the 
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X-ra y band, they have show n extreme spectral properties (see 
e.g.. ICos tamante et al.l 1200 lb and variability (e.g., Mkn 501, 
iPian et al . 1998). At very high energies (VHE, > 100 GeV), dou- 
bling timescales as short as a few minutes and flux variations 
of a factor of 10 in le s s than one hour h ave been observed 
jAharonian et al.l 12007 at lAlbert et al.ll2007l) . Their spectral en- 
ergy distribution (SED) is dominated by two broad peaks, lo- 
cated at UV-X-ray frequencies and in the GeV-TeV band. Their 
overall properties and behaviour have been most successfully 
- though not exclusively - explained as synchrotron and in- 
verse Compton (IC) emission from a population of relativis- 
tic electrons (see e . g.,|G hisellini et al. 1998; Spada et a l.ll200ll: 
Sikora & Madeiskil 12001: .Ghisellini et al . 2002; G uetta et al.l 
2004 and references therein), which upscatter their own 



self-produced synchrotron photons (synchrotron self-Compton, 
SSC: lKonigll [19811 iMaraschi et aP [l99l iBloom & Marsche^ 
Il996|) or external photons produced by different parts of th e 
jet (iGeorganopoulos & KazanasI I200I iGhisellini etal] |2005|) . 
Target photons ca n also be provided by the accretion disk and 
broad line region (Sikora et aIl ll994l;lGhisellini & Madaull 19961: 
Dermer & Schlickeiser 1993, BLR,) or by a dust y torus ( I Sikora! 
1994llWagner & Witzellll995l:ISikora et al.ll2002l) . In general, all 
of these seed photons can contribute significantly to the produc- 
tion of the observed SED, according to their energy density in 
the comoving frame. In HBL however, the lack of evidence of 
intense diffuse external fields (both directly from almost fea- 
tureless optical-UV spectra and from TeV 7-rays transparency 
arguments) has favoured the SSC model and external Compton 
process on photons from different parts of the jet as the most 
likely channels. 

Providing two handles on the one electron distribution 
responsible for both emissions, simultaneous observations in 
the X-ray and VHE bands represent both a powerful diag- 
nostic tool and a very str ingent testbed for the model itself 
(ICoppi & Aharonianll 19991) . especially during large flares when 
the emission from a single region is expected to dominate the 
SED. Alternatively, hadronic scenarios explain the 7-ray peak 
as being produced by ul tra high energy protons {E > lO^^eV; 
see e.g.lAharonia nll2000l: [Mucke & Protheroe 2001). 

Imaging atmospheric Cherenkov telescopes (lACT) provide 
a unique chance to study rapidly-variable sources at 7-ray en- 
ergies, thanks to their large collecting area. Multiwavelength 
campaigns performed on a few very bright sources (namely 
Mkn 501, Mkn 421 and lES 1959H-650) have shown that X- 
ray and VHE emission are generally highly correlated (e.g., 
Pianetall 119981: iDiannati-Atai etal] 119991; jAharonian et all 



1999blal; ISambruna et al.ll2000HKrawczvnski et al.ll2002L 120041: 



Blazeiowski et al .120051; iGiebelset al.l2007l:lFossati et al.i2008^ 
down to sub-hour timescales with no evidence of significant lags 
(IMaraschi et al.lll999l: [pbssati et al.ll2064l |2008|) . Moreover, the 
correlatio n seems to tighten when individual flares can be fully 
sampled (iFossati et al.ll2008l) . These results provide strong sup- 
port to the idea that both emissions during flares are produced 
by a single electron population ("one zone" SSC scenario). 

On the other hand, the same campaigns have also unveiled a 
more complex and puzzling behaviour, which represents a chal- 
lenge to the SSC scenario. Two main problems have recently 
emerged: 

1) the X-ray and VHE emissions do not always correlate 
(e.g., in Mkn 421; [Blazeiowski eTani2005l) . In particular, VHE 
flares seem to occur also without any visible X-ray counter- 
part (so-called "orphan" flares). The most striking example has 
been provided by lES 1959H-650 during the high state of 2002, 
when a strong (> 4 Crab) and rapid TeV flare (7 hours of 



doubling timescale) was not accompanied by det ectable varia- 
tions in the RXTE band ( iKrawczvnski et af]|2004l) . True orphan 
events are quite difficult to explain with a standard SSC scenario. 
However, the generally sparse sampling does not allow the ex- 
clusion of lagged counterparts dBlazeiowski et al.l2005b or coun- 
terpar ts emerging in a different energy band ( Krawczvnski et al.l 

[lool . 

2) Mkn 421 exhibited a quadratic relation between VHE and 
X-ray flux variations during both the rising and decaying phases 
of a flare ( iFossati et al.ll2004 l2008l) . This is not expected if the 
source is in the Klein-Nishina (KN) regime. 

A blazar is generally said to be "in the Klein-Nishina 
regime" when the observed 7-ray emission is produced by TeV 
electrons which do not upscatter their own self-produced syn- 
chrotron photons - since inhibited by the smaller cross-section 
of the KN regime - but upscatter in the Thomson regime lower- 
energy photons produ ced by lower-energy electron^ (see e.g., 
iTavecchio etanil998l) . This condition changes the mapping of 
the synchrotron and Compton components, so that the two peaks 
are not produced by electrons of the same energy. In this situ- 
ation, the VHE emission tends to track the X-ray synchrotron 
variations only linearly instead of quadratically, although a flare 
that is achromatically extended over a sufficiently wide energy 
range can still yield a quadratic increase (iFossati et al. 20Qa). 
However, the energy dependence of both synchrotron and IC 
cooling (cx 7^) prevents this relation during the decaying phase: 
since higher-energy electrons cool faster than lower-energy elec- 
trons, they see a roughly constant seed-photon energy density, 
resulting in a mostly linear decrease. A quadratic decrease could 
be achieved by imposing the strict Thomson condition, but that 
seems to requ ire extrem ely large beaming f actors for Mkn 421 
dKatarzvhski e t al. 20051; IFossati et al.ll2008l) . 

To investigate these issues, a multiwavelength study of sin- 
gle flares is essential. Although many efforts have been made to 
achieve a good sampling, so far the short variability timescales 
have been difficult to study, because of the lack of sensitivity of 
the past-generation of lACTs. However, these are extremely in- 
teresting timescales: the results of HEGRA on Mkn 421, for ex- 
ample, have already revealed an ent ire "zoo" of intra-night fl ares 
with different rise and decay times dAharonian et alj[2002bh . in- 
dicative of a complex inte rplay between acceleration/ injection 
and cooling processes (e.g.. lKirk & Mastichiadisiri999t) 

When the lACT array HESS became operational, a project 
was therefore developed with specific ToO proposals, to investi- 
gate the fast variability timescales with a whole night (6-8 hours) 
of continuous, simultaneous observations during a bright 7-ray 
state. To achieve this aim, Chandra was chosen because it is 
the only X-ray satellite capable of a full coverage of the entire 
HESS visibility window during most of the year, and without 
the interruptions on sub-hour timescales which are typical of 
low-orbit satellites. The efforts paid off in July 2006, when the 
HBL PKS 2155-304 (z=0.116) became highly active at VHE, 
with a flux level a factor of ~ 10 higher than its typical quiescent 
flux of -4 X IQ-ii cm-2s-i above 200 GeV. PKS 2155-304 is 
one of the brightest and most studied BL Lacs in the Southern 
Hemisphere, at every wavelength, and it can be detected by 
HESS on almost a nightly basis since 2002 (see lAharonian et aU 
2005a b, and references therein). 



' Note that this case is different from the condition where the ob- 
served 7-rays are indeed produced by IC-scatterings occurring in the 
Klein-Nishina regime, or w hen the cooling itself is determined by 
Klein-Nishina losses, see e.g., iModerski etaP l l2005l) . 
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In the first days of activity, ToO observations were also trig- 
gered for other X-ray satellites such as RXTE and Swift, to sam- 
ple the source behaviour over several days and weeks. Then, in 
the early hours of July 28 2006, a giant outburst occurred at 
VHE (^lOOx above the quiescent level), with a peak flux of 
15 Crab above 200 GeV (corresponding to 9.9 Crab above 400 
GeV ) and repeated flares with doubling timescales of few min- 
utes dAhar onian et al. 2007a). Unfortunately, this dramatic out- 
burst occurred too early with respect to the already-triggered but 
not-yet-started X-ray observations (which were acquired from 
the night after). Therefore, the event of most exceptional vari- 
ability remained without multiwavelength coverage. 

Two days later however, on the night of July 29-30, the 
source underwent a second major VHE flare, this time in co- 
incidence with our planned C/iant/ra-HESS ToO campaign, and 
with the further coverage in the optical band provided by the 
Bronberg Observatory in South Africa. Snapshot observations 
of few ks were also taken with RXTE and Swift. This second out- 
burst has reached even higher fluxes than the first one (~ 11 Crab 
above 400 GeV). As a result, in this single night the most dense 
and sensitive X-ray /TeV campaign to date was obtained, during 
one of the brightest states ever observed from an HBL at VHE. 

This paper focuses on the multiwavelength results of this ex- 
ceptional night, presenting new optical. X-ray and VHE data. 
The results of the whole VHE activity of PKS 2155-304 be- 
tween July and October 2006 will be presented in forthcoming 
papers, together with the overall multiwavelength coverage. 

Since this is a very rich and complex dataset, the data have 
been divided into several subsets of different time windows, to 
highlight specific aspects (e.g., different VHE thresholds, inte- 
gration times, or X-ray coverage). These subsets and their ratio- 
nale will be introduced in the relevant Sects., but a summary list 
with corresponding time windows is given in Table [T] for ref- 
erence. Throughout the paper, the following cosmological pa- 
rameters are used: Hq = 70 km s^^ Mpc^^, f^M = 0.3, and 
ft\ ~ 0.7. Conforming to the convention adopted in all pre- 
vious HESS publications, unless otherwise indicated, all errors 
are given at the 1 <t confidence level for one parameter of interest 
(Ax^ = 1.0). For simplicity, in the text the values of MJD are 
given as MJD*= MJD-53900. 

2. Observations and data reduction 

2.1. HESS 

HESS is an array of four Imaging Atmospheric Cherenkov 
Telescopes located in the Khomas Highlands of Namibia (23°S, 
15°E, 1800 m a.s.l.). Each telescope has a surface area of 107 
m^ and a total field of view of 5°. A camera consisting of 960 
photo-multipliers is located at the focal length of 13 m. Each 
camera images the dim Cherenkov flashes from air-showers of 
VHE 7-rays in the atmosphere, collected by its mi rrors (for more 
details about the layout of the telescopes, see iBernlohr et al.l 
(l200l ). The camera images are ca librated following the pre- 
scriptions in lAharonian et al.l (l2004h . The stereoscopic view of 
the air showers allows the reconstruction of the direction of the 
primary 7-ray with an accuracy of w0.1° following method 1 of 
iHofmann et all d 19991) . 

The recorded signal in the field of view is dominated by the 
constant background from hadronic cosmic rays entering the at- 
mosphere. Most of the hadronic background can be identified 
from by the shape of the shower images and the arrival direction 
of the recorded showers. The remaining hadronic background 
can be statistically removed by estimating it from sky regions 



with no 7-ray signal. In the analysis shown here, loose cuts 
and the reflected background method were applied for back- 
ground substraction (jXharonian et al. 2006bJ. Light curve and 
spectra were derived following the standard HESS analysis also 
described in this reference. 

HESS observed PKS 2155-304 throughout the night of July 
29-30. A total of 15 runs (each 28 min long) were taken, all 
passing the standard quality criteria specified in lAharonian et al.l 
( l2006bh . The total Ufetime after dead-time corrections is 6.58 
hours. A 7-ray excess of 32612 events was dete cted with 
a sign ificance of 254 a following Equation 17 in iLi & Mai 
(Il983h . The excess is point-like, taking into account the point 
spread function of HESS, with a best fit position of Q!2ooo = 
2l''58'"52.6" ± OA%^f ± 1.3^^^,, (Saooo = -30°13'29.8" ± 
1.9'^tat ± 20'4,, consistent with the position of PKS 2155- 

304 (aannn = 2l''58'"52.065P, (^2000 = -30°13'32.118"; 
iMa et alj|l998l) . Because of the long duration of the observa- 
tions, the zenith angle of the source varied strongly during the 
night, going from 53° at the beginning, to 8° at the middle, 
and back to 50° at the end of the night. Observations at larger 
zenith angles imply a higher energy threshold of the analysis. 
For the observations discussed here the energy threshold varies 
between 200 and 7 00 GeV for the applied energy reconstruc- 
tion (.Aharonian et al . 2006il). This results in a tradeoff between 
energy and time coverage for the analysis, as we discuss later. 

HESS has a systematic uncertainty in the normalization of its 
energy scale of wl5%. The main source of this systematic error 
are uncertainties in t he atmospheric conditio ns (for a more de- 
tailed discussion see lAharonian et aT]|2006bl) . During the night 
the atmospheric conditions were stable. This can be verified to 
timescales shorter than one minute in the overall trigger rate and 
background rates in regions off the source. The differential en- 
ergy spectrum ^{E) of PKS 2155-304 at VHE energies is gen- 
erally steep, with a photon index of about 3.4 in a power-la w 
model <S>{E) = dN/dE = <i>oE-^ ( lAharonian et al.ll2005allR 
l2007al) . The systematic error in the energy scale of the detec- 
tor therefore transforms into an error of w40% in the overall 
flux normalization. The systematic error in the slope of dif- 
ferential energy spectr a is wO.l for the photon index F (see 
lAharonian e"tani2006bl) . 

2.2. Chandra 

PKS 2155-304 was observed with Chandra dWeisskopf et al.l 
I2OOOI) for a total duration of 30 ks with the Low Energy 
Transmission Grating (LETG) spectrometer coupled with the 
ACIS-S detector (Obsid 6874; set-up with 1/8 subarray and 6 ac- 
tive chips, for a 0.7 s frametime). Because of difficulties in plac- 
ing our constrained ToO observation within the Chandra sched- 
ule, the observation started later than requested, missing the first 
1.7 hours of the HESS window. Unfortunately, the main flare oc- 
curred in the first few hours of the HESS window. As a result, 
the rising part of the main VHE flare has no X-ray coverage. 

Data reduction was performed according to the standard 
CXC procedures, using the CIAO software version 3.4 with the 
corresponding Calibration Database CALDB version 3.3.0 and 
HEASOFT v6.3.2. Event files on timescales as short as 2 min- 
utes were obtained using dmcopy, which propagates all dead- 
time corrections correctly. Grating spectra were then extracted 
with tgextract and their ancillary files were generated with 
fullgarf for each arm; then added together to obtain the first- 
order spectrum. The scientific analysis was completed mainly on 
the first-order spectrum, because of its higher S/N and photon 
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Statistics. A check was performed that the centroid of the source 
obtained from the zero-order image was indeed coincident with 
the coordinates of the source on the detector The response ma- 
trix was produced using mkgrmf applied to the entire observa- 
tion, since no difference was found from files created in different 
epochs during the night. The analysis was optimized and per- 
formed only on the continuum properties: the study of the total 
grating spectrum at its highest resolution is beyond the scope of 
this paper and is left to future publications. The background and 
source spectra in each time-bin were then extracted with the tool 
tg_bkg for use in XSPEC. 

The hundreds of spectra (one for each time bin, down to 2- 
minute bins) have been routinely fitted in XSPEC version 11 .3.2, 
using a source model plus photoelectric absorption (wabs), with 
the equivalent hydrogen col umn density fixed at Galact ic values 
(A^H = 1-69 X 10^" cm~^; iDickev & Lockmanlll990h . This is 
also the A^h value obtained from the best fit to the total exposure 
with free absorption, to within 1 sigma. The integrated flux and 
its error were calculated from the spectral fit using the specific 
error routine in XSPEC. The error in the unabsorbed flux was 
then obtained from the percentage error of the absorbed fit. The 
results of this procedure were later checked to be fully consistent 
with the values from the specific Tel routine f luxerror . tcl 
recently provided by HEASARC with XSPEC vl2, for calculat- 
ing the error in the flux from single components of the modeQ. 
In the following, all X-ray fluxes are quoted as unabsorbed val- 
ues. The time analysis was also performed using the direct count 
rate and its error in the energy band of interest, for each time- 
bin, obtaining fully consistent results. The average count rate 
observed from PKS 2155-304 in the LETG is 8 cts/s, ranging 
between 12 and 6 cts/s. These count rates allow the spectra to be 
extracted down to 2-4 minute timescales with typically 1000- 
2000 counts each. For the observed flux, the grating spectra do 
not suffer from pile-up problems: the total fraction of piled-up 
events estimated at the peak of the effective area (1-2 keV) and 
source peak flux is less than 5%. 

One of the calibration issues with the ACIS instrument is 
the excess absorption seen below 1 keV due to the build-up 
of contaminants on the optical blocking filter These contami- 
nants (thought to be carbon compounds) cause a significant ab- 
sorption in the 0.3-0.4 keV range, which is taken into account 
by the calibration but which also severely reduces the count 
rate in that range, yielding very few counts during short expo- 
sures. Therefore, for the spectra extracted on 2 and 4 minute 
timescales, the interval 0.3-0.4 keV was excluded from the anal- 
ysis. At high energies, data were included in the fit up to the en- 
ergy where positive net source counts were present. As a result, 
short-timescale spectra were fitted in the range 0.23-0.3 and 0.4- 
6 keV, while longer-exposure spectra could be fitted from 0.23 to 
8-9 keV. The spectra were generally rebinned to have more than 
30 counts per bin, using different schemes. A fixed coarse bin- 
ning was used for all the spectra on short timescales (< 7 min- 
utes). Various checks have shown that, within the uncertainties, 
the obtained results are independent of the adopted rebinning. 

The analysis of the total Chandra exposure shows evidence 
of a slight residual excess absorption in the 0.3-0.4 keV range, 
not yet fully accounted for by the calibration. The effect is small, 
and does not affect the fit values significantly. Nevertheless, 
it was taken into account by simply adding an edge model at 
the carbon K edge, 0.31 keV, with a fixed value of r = 0.4. 
These are the best-fit parameters obtained from the fit of the to- 
tal Chandra spectrum (see Sect. 14. 2t . 



2.3. RossiXTE 

As part of the ToO campaign for a daily monitoring, RossiXTE 
dJahoda et alJll996l) pointed PKS 2155-304 twice during the 
night of July 29-30, for a total exposure of ^2.6 ks (44 min- 
utes). Fig.[2]shows the epochs of the RXTE windows with respect 
to the overall HESS and Chandra light curves, together with the 
time-windows of the Swift obs e rvatio n. The latter is analyzed 
and discussed in iFoschini et al.l (|2007|) . Since the Swift pass- 
band mainly overlays the Chandra band, we focus only on the 
RXTE spectrum, to extend the X-ray spectrum in the hard band. 
The data reduction and analysis were performed using FTOOLS 
v6.3.2 with the standard procedures and filtering criteria recom- 
mended by the RXTE Guest Observer Facility after September 
2007 0. For a more accurate spectral determination, only the 
PCU2 data were considered. The average net count rate from 
the source is measured at 6 cts/s/pcu, in the 3-20 keV band. 

The RXTE spectrum was then fitted in XSPEC together with 
the Chandra spectrum extracted from the same time window, ob- 
taining a spectral measurement over two decades in energy (0.2- 
20 keV). Without adjustments, the two spectra have very similar 
normalizations, indicative of a very good inter-calibration be- 
tween the two instruments. To obtain the most accurate spectral 
determination, the RXTE /Chandra normalization was fixed at 
the value measured by fitting the same power-law model in the 
overlapping energy range (3-7 keV), namely 1 .08 (see Sect. l4.2b . 

2.4. Optical data 

Optical observations in the V filter were performed using the 
32 cm Schmidt-Cassegrain telescope at Bron berg Observatory, 
Pretoria, South Africa (see llmada et al.ll2008l for details about 
observations of variable sources with this telescope). After de- 
biasing and flat fielding, data were analyzed using relative aper- 
ture photometry with a K-type star of similar magnitude (12.6). 
Frames were taken every 30 s and then smoothed over 6 succes- 
sive data points to calculate the mean value in each bin. Each 
point of the optical light curve therefore has a time duration of 
-180s. 

The comparison with the reference star shows that the rise in 
optical intensity is highly significant. The K-star light curve is 
constant with an intrinsic scatter in the datapoints - using the 
same smoothing procedure - of <0.02 mag. The errors were 
determined from the variance of each 6 successive data points 
for both the source and reference star. No time variability in the 
intensity of the reference star is seen. The optical fluxes were 
corrected for Galactic extinction with Av=0.071 mag. 



3. Temporal analysis 

3.1. TeV light curves 

The measured VHE light curve in two-minute time bins is shown 
in Fig.[ri divided into five different energy bands. The time cov- 
erage increases with threshold energy, due to the aforementioned 
zenith-angle effect. A full coverage of the observation is thus 
reached for an energy threshold of 700 GeV, while light curves 
(and spectra) down to —300 GeV are obtained only for the cen- 
tral five hours of the observing window (see Table [T]). The dif- 
ferent time windows have therefore been labelled according to 



^ See |http://xspec.gsfc.nasa.gov/docs/xanadu/xspec/fluxerror.htnil 



^ See |http://www.universe.nasa.gov/xrays/programs/rxte/pca/doc/bkg/bkg-2007^ 
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Fig. 1. VHE fluxes integrated in different energy bands, as a function of time, in time bins of two minutes. The time windows 
corresponding to the different energy thresholds are given in Table [T] (labelled accordingly, from top to bottom: T200 to T700). 
The dotted lines mark the positions of rapidly varying events on-top of the main flare (see text). 



the respective energy thresholds (from T2 to T7 0), for ref- Nebula flux dAharonian et al.ll2006bl) above the same threshold, 

erence. This peak flux is about 20% higher than the peak flux measured 

The source underwent a major flare in the first hours of ob- during the night of July 27-28, above the same threshold (~9 

servation, reaching a peak flux of -9.9 x lO^^" cm'^s"! (>400 Crab, -8 x lO^^" cm^^g-i >400 GeV). The fluxes measured 

GeV) at MJDM5.90, corresponding to -11 times the Crab in these two nights (July 27-28 and 29-30) are the highest ever 
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Table 1. Summary of the subsets of VHE data used in this paper (MJD*=MJD-53 900). 



Label" 




MTD^start 


MtD*enfl 


ration 

[hrs] 


[GeV] 




Sect.*" 


notes 


T200 




45.934618 


46.107844 


4.16 


200 


liJ 


§3.1 




T300 




45.913536 


46.129121 


5.17 


300 


□ 


§3.1 




T3 0- 




J . "Z. JZ. JZ. 


46 179171 


4.94 


300 


l6l 


83 3 


T3 with X~rfly covGrsgG 


T400 




45.892483 


46.150125 


6.18 


400 




§3.1 




T500 




45.871450 


46.171258 


7.20 


500 


□ 


§3.1 




T700 




45.850393 


46.171258 


7.70 


700 


□ 


§3.1 




T4 0- 


-Peak 


45.896643 


45.920474 


0.57 


400 


|8J 


§4.1 


Spectrum at the peak of VHE flare 


T300- 


-High 


45.913530 


45.970312 


1.36 


300 


m 


§4.1 


Average high-state spectrum inside T300. 


T300- 


-Low 


46.013252 


46.129166 


2.78 


300 


m 


§4.1 


Average low-state spectium inside T300. 


T300- 


-Xmax 


45.922130 


45.944699 


0.54 


300 


m 


§4.2 


Highest simultaneous X-iay/VHE state. 


T400- 


-Xmin 


46.100000 


46.150000 


1.20 


400 


m 


§4.2 


Lowest simultaneous X-ray/VHE state. 


T300- 


-RXTE 


46.016846 


46.033328 
















46.084624 


46.098698 


0.73 


300 


m 


§4.1 


Sum of the 2 intervals with RXTE coverage. 



" The T-label corresponds to the VHE energy threshold of that dataset. 

* Figure where the time intervals are indicated and the Sect, where they are first introduced. 



observed at VHE from PKS 2155-304. The total amplitude of 
flux variation during this night was about a factor of ~20, above 
400 GeV (T4 covers both the highest and lowest flux epochs), 
similar to the flux variation observed in the night of July 27-28 
(^23x). 

The main flare seems to occur with similar rise and de- 
cay timescales, of the order of 1 hour (half-to-maximum am- 
pli tudes, measured using a "generalized Gaussian" function as 
in lAharonian et al.l l2007ah . After the peak, the VHE flux de- 
creased overall during the night reaching its minimum around 
MJD*46.12, but with two other smaller-amplitude flares super- 
imposed: a short burst around MJD*45.96 of duration ~20 min- 
utes, and a longer flare or plateau between MJD*46.0 and 46.1, 
with a duration of 2-3 hours. 

In addition, two further sub-flares are evident in all cov- 
ered energy bands, around MJD*45.885 and MJD*45.920 (dot- 
ted lines in Fig.IIli- These structures have a duration of ~10 min, 
similar to the flares of the night of July 27-28. Although there 
are hints of even shorter variability (few minutes), the signifi- 
cance is limited. 

3.2. Comparison with X-ray and optical liglit curves 

The combined VHE, X-ray and optical light curves are shown 
in Fig. |2] and [3] Significant flaring activity is visible in all three 
bands, but with different amplitudes. To emphasize the specific 
variability patterns, the vertical scales in Fig. |2] were adjusted 
differently for each band. Fig. |3] shows instead the light curves 
on the same flux scale, but with a i^F,^ representation. They cor- 
respond to slices of the SED at the three energies of 0.3 TeV, 
0.3 keV, and 2.25 eV (i.e. 5500A). In this way, it is possible to 
highlight the overall changes and time evolution SED-wise. 

The 0.3 TeV fluxes were calculated from the integrated 
>300 GeV light curve (T300 window) using the average 
power-law spectrum measured in the respective epochs (namely, 
the T300-High spectral index in the high state, and the 
T30 0-LOW index elsewhere, see Sect. 4.1). The same proce- 
dure was used to calculate the 0.3 TeV fluxes from the >500 
GeV light curve (T5 00 window), during the epoch not covered 
by T3 (empty circles in Fig.[3]i. A comparison with the results 
of the >300 GeV light curve in the overlapping window shows 



that the extrapolation from 500 GeV does not introduce differ- 
ences of more than ~2%. The 0.3 TeV fluxes were then cor- 
rected for the absorption effect caused by 7-7 interactions with 
the Extragalactic Background Light (EBL), using the model by 
Franceschini et al. (2008) (discussed in Sect. 4.1.1). This corre- 
sponds to a low density of the EBL, close to the lower limits 
obtained by galaxy counts. The plotted fluxes therefore can be 
considered as lower limits to the intrinsic VHE emission of the 
source. The X-ray fluxes at 0.3 keV were calculated with the 
same procedure, using the power-law spectrum fitted in each of 
the individual bins. Anticipating the result that both the VHE and 
X-ray spectra are steep (r>2; see Sect. 4), the plotted i/F^ fluxes 
- at the low-energy end of the respective passbands - provide an 
estimate of the emission closer to the respective SED peaks than 
the i/Fi, fluxes in the hard band. 

Several remarkable features can be noted. The first is the 
huge difference in amplitude between the variations in the three 
energy bands. In few hours, the VHE flux changed by more than 
an order of magnitude, whereas the X-ray flux varied by only a 
factor ^2 and the optical flux by less than 15% (the contribution 
of the host galaxy is negligible, see Sect. 6.3). The source thus 
shows a dramatic increase in variability with photon energy. 

Secondly, the VHE emission dominates the energy output 
by far in the three bands. In Fig. [3] the comparison between the 
0.3 TeV and 0.3 keV fluxes shows the evolution of the Compton 
dominance of the source, i.e. the ratio of the IC 7-ray luminosity 
to the synchrotron power (Lc/Ls). The 7-ray luminosity dom- 
inates the synchrotron luminosity by a factor ^^8 close to the 
flare maximum, evolving rapidly towards comparable levels at 
the end of the night. This is the first time that such high and 
rapidly variable Compton dominance is observed in an HBL, ir- 
respective of the choice of EBL density. The swiftness of the 
changes in the 7-ray emission and SED properties also under- 
lines the danger in modeling X-ray and VHE data taken even 
only a few hours apart, during such events. As shown in Fig.|2] 
both RXTE and Swift observations occurred when the Compton 
dominance had already decreased significantly. 

Thirdly, despite this difference in amplitude, the X-ray and 
VHE light curves are strongly correlated, with the X-ray emis- 
sion following closely the same pattern traced by the VHE light 
curve (see next Sect.). The optical light curve, instead, does not 
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Fig. 2. Overall light curves of PKS 2 1 55-304 in the night of July 
29-30 2006, as seen by HESS (T7 00, upper panel), Chandra 
(lower panel, blue circles), and the Bronberg Observatory (op- 
tical V band, red squares). Time bins of 4 minutes (3 for the 
V band). The segments on the upper x-axis also show the two 
intervals corresponding to the RXTE exposure (R label), and 
the two intervals of the Swift pointing (S label) reported in 
iFoschini et al.l (l2007h . The vertical scales differ in each panel, 
and have been adjusted to highlight the specific variability pat- 
terns. Lower panel: the left axis gives the integrated 0.5-5 keV 
flux, the right axis gives the V-band i/Fi, flux at the effective fre- 
quency 5500 A; both are in units of erg cm~^ s~^. The vertical 
line marks a visual reference time for the start of both the op- 
tical and VHE flares. The shaded bands mark the time interval 
where the T300-Xmax and T400-Xmin dataset are extracted 
(highest and lowest X-ray/VHE state; see Table 1). 



correlate on short timescales with the other two bands. However, 
there is a rise of ~15% in flux, which appears to begin at the 
same time as the main VHE flare. A conservative estimate of the 
chance probability of coincidence -considering only the data of 
this night- is of the order of a few percent. We discuss in Sect. 
8.2 the possible implications if this simultaneity is genuine. 

3.3. Inter-band time lags 

The degree of correlation and possible time lags between dif- 
ferent light curves have been quantified by means of cross- 
correlation functions. The correlation analysis was performed 
between X-ray and 7-ray light curves and between hard and soft 
energy bands within each passband. As mai n tool it was used 
the dis crete correlation function (DCF) from l Edelson & KrolikI 
(Il988h . The DCF is especially suited to unevenly spaced data, 
such as the VHE light curves, which have gaps of a few min- 
utes every 28 minutes between the stop and start of consecutive 
runs. The time lags between light curves are determined to be the 
maximum of a Gaussian plus linear function fitted to the central 
peak of the DCF. This procedure is robust agai nst spurious peaks 
at zero time-tag caused by correlated errors (lEdelson & Krolilj 



S -9 




- 46.5 :r 



45.8 



45.9 



46 46.1 
MJD-53900 



46.2 



Fig. 3. Light curves of the uF^ flux at 0.3 TeV, 0.3 keV, and 
55OOA in 4-minute bins, plotted on the same flux scale. The 
right axis reports the corresponding luminosity scale. Note the 
remarkable difference of the amplitudes and the dramatic evolu- 
tion in the VHE/X-ray flux ratio. HESS data (black): filled cir- 
cles are fluxes calculated from the T 3 light curve (integrated 
above threshold); empty circles from the T5 light curve (see 
text). The 0.3 TeV fluxes are corrected for intergalactic 7-7 ab- 
sorption with a low-density EBL model, while both X-ray and 
optical fluxes are corrected for Galactic extinction (Av=0.071). 



Il988h . The error in the measured time lag is determined by sim- 
ulations. Ten thousand light curves were generated by varying 
each measured point within its errors. The entire correlation pro- 
cedure was repeated for each of these simulations, resulting in 
a cross correlation peak distribution (CCPD). The RMS of the 
CCPD then provides an estimate of the statistical error in the 
measured time lag (Maoz & Netzer 1989; Peterson et al. 1998j). 
The time binning of the light curves and the DCF introduces an 
additional systematic error. The latter was estimated to be 30 s 
by injecting various artificial time shifts into the original VHE 
photon list, smaller than the duration of the time bins, and by 
measuring the relative shift of the CCPD. 

At VHE, the analysis was performed on the simultaneous 
light curves between 300-700 GeV and above 700 GeV in two- 
minute time bins (Fig.[T]). This choice yields a good compromise 
between event statistics, time coverage, and a maximum energy 
difference between the bands. The resulting time lag between 
the higher and lower energy band is (28 ± 30stat ± 30sys) s (see 
Fig.lHi. This time lag does not dif fer significantly from ze ro, as 
for the flare of two nights before ( lAharonian et al.ll2008bh . and 
we derive a 95% confidence upper limit of 129 s. This value 
was calculated by assuming a Gaussian probability distribution 
around the measured time lag. The width of the distribution was 
set to be the linear sum of the statistical and systematic error, 
to be conservative. Afterwards symmetric intervals around zero 
were integrated, until a 95% containment was achieved. 

In the X-ray band, an analogous procedure was applied. The 
total Chandra light curve was divided into a soft (0.2-1.0 keV) 
and hard (2.0-6.0 keV) band. Because of the larger errors, in 
this case 4-minute time bins were used. The measured time lag 
is (-82 ± 202stat ± 30sys ) s. This value again does not differ 
significantly from zero, resulting in a 95% upper limit of 482 s. 
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Fig. 4. Cross-correlation analysis of the hard versus soft bands at 
VHE (upper panel) and at X-ray energies (lower panel). Upper 
panels: (a) DCF of the 300-700 GeV and the >700 VHE band. 
The line around the peak shows the best fit Gaussian plus lin- 
ear function, with a maximum at 28 s. (b) Corresponding cross- 
correlation peak distribution (CCPD) of 10000 simulated light 
curves. The RMS of the distributions is 30 s. The dotted line 
marks the position of the maximum in (a). Lower panels: (c) 
DCF of the 0.2-1.0 keV and 2.0-6.0 keV X-ray band. The line 
around the peak shows the best fit Gaussian plus linear function, 
with a maximum at -82 s. (d) Corresponding CCPD of 10000 
simulated light curves. The RMS of the distributions is 202 s. 
The dotted line marks the position of the maximum in (c). 
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Fig. 5. Cross-correlation between the X-ray and VHE light 
curves: (a) DCF of the >300 GeV light curve and the LETG 0.2- 
6 keV band. The blue line shows the best fit Gaussian plus linear 
function, with a maximum at -10 s. (b) Corresponding cross- 
correlation peak distribution of 10000 simulated light curves. 
The RMS of the distributions is 76 s. The dotted line marks the 
position of the maximum in (a). 



The cross-correlation analysis between the X-ray and 7-ray 
emission was performed on the simultaneous light curves with 
two-minute time bins shown in Fig. |6] (in the T30 0-X time 
window). The resulting cross-correlation is shown in Fig. |5] 
The two light curves overall are highly correlated, with a max- 
imal correlation of DCFmax ~ 0.9, and no significant lag is 
found. The time lag of the X-rays with respect to the 7-ray is 
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Fig. 6. HESS and Chandra light curves in the simultaneous time 
windows corresponding to a 300 GeV threshold (T300-X in 
Table 1). Two-minute time bins. Upper panel: integral fluxes 
above 300 GeV. Lower panel: Ist-order LETG count rate in the 
0.2-6 keV band. 

(-10 ± 76stat ± 30sys) s, yielding a 95% confidence upper limit 
to a time lag of 208 s. To test whether this result was caused by 
an averaging of lags with opposite signs, the correlation analy- 
sis was also performed on sub-intervals, namely in the interval 
around the first smaU flare at MJD*45.96 (MJD*45.94-46.0) and 
after MJD*46.0. The two emissions are again highly correlated 
(DCFmax ~ 0.9), with no evidence of time lags. The strong cor- 
relation is determined not merely by the overall decaying trend 
of both light curves, but also by their specific patterns: a DCF 
max value of w 0.7 — 0.8 is still obtained after whitening the 
light curves by removal of either a linear or quadratic trend. On 
the shortest timescales (<4-8 minutes), however, the VHE light 
curve shows few small flares apparently not mirrored in the X- 
ray band (see e.g. MJD*45.925 and 46.060 in Fig.|6ll. Although 
this might indicate a more complex correlation on the fastest 
timescales, at present no firm conclusions can be drawn, since 
the significance of these structures is low. 

The correlation analysis was also performed using the task 
CROSSCOR of the Xronos 5.21 package, which measures the 
correlation function (CCF) with a direct Fourier algorithm. This 
algorithm requires a continuous light curve without interrup- 
tions, therefore the few gaps were filled with the running mean 
value calculated over the 8 closest bins (e.g., iRavasio et al.l 
,2004t ). The results are fully consistent with the DCF analysis, in- 
dicating that the small gaps in the VHE light curves - the X-ray 
light curve is continuous - do not introduce significant distor- 
tions for such well sampled data. 

The cross-correlation analysis between VHE and X-ray light 
curves was limited to the strictly simultaneous window, to avoid 
artifacts in the lag determination due to the different timespans 
and the light curves characteristics. Because both light curves 
have each one main flaring feature, the cross-correlation per- 
formed on different intervals tends simply to match the max- 
ima of the two emissions in those intervals, irrespective of 
the smaller amplitude patterns. This would yield an artificial, 
window-dependent "lag" with typically lower correlation values 
(as is the case here, with a timespan of ^^3200 s between the 
maxima of the overall VHE and X-ray light curves and lower 
DCF/CCF values -0.7). 
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Fig. 7. Left panel: integral flux >300 GeV and photon index as a function of time (T300 dataset). Horizontal error bars show the 
time interval of each bin, going from 7 to 14 minutes before and after MJD*46.0. The shaded zones mark the two time-intervals 
corresponding to the average high and low-state spectra fitted in Table|2](T3 0-High and T3 0-Low, respectively). Right panel: 
photon index as a function of the integral flux. 
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Fig. 8. Left panel: integral flux >500 GeV and photon index as a function of time (T5 00 dataset). Horizontal error bars show the 
time interval of each bin, going from 14 min. at the beginning to 28 min. towards the end of the night. The shaded area shows the 
time window where the 7-ray peak spectrum (T4 00-Peak dataset in Tables [T|2] | has been extracted from. Right panel: photon 
index as a function of the integral flux. 



4. Time-resolved spectral analysis 

4.1. VHE spectra 

A search for spectral variations in the VHE data was performed 
by fitting a power-law spectrum to a fixed energy range in fixed 
time bins. The unprecedented statistics of this dataset allow the 
sampling in 7 to 1 4-minute bins in the T 3 time window (Table 
1). On these short integration times, the power-law function 
gives a statistically good description of the data. The results are 
shown in Fig.|2] The spectrum generally hardens with increasing 
flux. The fit to a constant photon index results in a probability 
of only 1.6%. These spectra are also used for comparison with 



the X-ray spectra extracted in exactly the same time bins, and 
discussed in Sect.|5] 

A study of the spectral variations was also performed on the 
T5 00 dataset, which allows the sampling of a wider time span 
and in particular of both the rise and decay phases of the main 
7-ray flare. However, the lower number statistics requires longer 
integration times, yielding a lower time resolution. The result 
is shown in Fig. |8] where spectra were extracted in 14 and 28 
minutes bins. The spectral variations follow the same pattern as 
for the T3 spectra, both in time and in the flux-index relation. 
No significant spectral changes are observed between the rising 
and decaying part of the flare, with the possible exception of a 
hardening event that precedes the peak of the 7-ray emission by 
«28 minutes. 
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Fig.9. Selected VHE spectra: average spectra for the T300-Low (left) and T300-High (middle) datasets (see Fig. |7]i, and the 
average spectrum above 200 GeV (T2 0, right). Open symbols corresponds to the spectra corrected for EBL absorption as described 
in Sect. 4.1.2. 



Table 2. Spectral fit of the measured VHE spectra extracted in different epochs (see Table 1). 



PL fits 


$0 




r 






^0.3-3 ToV 


Xr (d.o.f.) 






cm^^s'^TeV" 


1 








erg cm~^ s~^ 






T400-Peak 


2.20±0.06 X 10' 


-10 


3.46±0.04 






1.35 xlQ-"" 


3.66(11) 


3 xlO"^ 


T300-High 


1.36±0.03 xlQ- 


-10 


3.36±0.03 






7.88 xlQ-^" 


4.3(12) 


7 xl0~^ 


T300-LOW 


4.89±0.15 xlO' 


-11 


3.51±0.03 






3.09 xlQ-^" 


1.36(12) 


0.18 


T200 


7.46±0.12 xlO" 


-11 


3.25±0.01 






4.06 xlO""' 


16(16) 





T300-RXTE 


4.78±0.30 XlO" 


-11 


3.53±0.07 






3.07 xlQ-^" 


1.4(11) 


0.17 


PL exp. cut. 


$0 




r 


£cut(TeV) 




-fo.3-3 ToV 


Xr (d.0.f.) 


P.^ 


T400-Peak 


5.96±1.19 xlO" 


-10 


2.60±0.17 


1.19±0.25 




1.28 XlO""" 


0.38 (10) 


0.96 


TSOO-High 


2.72±0.36 XlO" 


-10 


2.86±0.09 


1.66±0.33 




7.95 xlO"^" 


1.41 (11) 


0.16 


TSOO-Low 


7.08±1.25 XlO- 


-11 


3.26±0.12 


2.92±1.42 




3.08 xlQ-i" 


0.98 (11) 


0.47 


T200 


2.12±0.23 XlO- 


-10 


2.65±0.06 


l.OliO.ll 




4.25 xlO"^" 


4.41 (15) 


2 xlO"* 


Log-parabolic 


$0 




r 


b 




-Fo.3-3 TcV 


Xr (d.0.f.) 


P.^ 


T400-Peak 


2.55±0.09 XlO- 


--10 


3.54±0.06 


1.05±0.20 




1.25 xlO""^ 


0.37 (10) 


0.96 


TSOO-High 


1.46±0.04 xlO- 


-10 


3.53±0.05 


0.62±0.11 




7.94 xlO"^" 


1.35 (11) 


0.19 


T300-LOW 


4.98±0.16 xlO- 


-11 


3.66±0.07 


0.41±0.16 




3.08 xlO"^" 


0.76(11) 


0.68 


T200 


7.51±0.16 xlO- 


-11 


3.69±0.05 


0.78±0.07 




4.29 xlO"^" 


2.39(15) 


0.002 


Broken PL 


$0 




Ti 


Ta 


-Bbroak(TeV) 


-fo.3-3 ToV 


xl (d.o.f.) 


P.^ 


T200 


1.46±0.10 xlO- 


-10 


2.73±0.05 


3.60±0.04 


0.42±0.02 


4.31 xlO"^" 


1.44 (14) 


0.12 



The study of the spectral shape in more detail requires 
higher event statistics. To achieve this, the dataset was divided 
into similar spectral states, namely a high and a low flux state 
(T300-High and T300-Low, respectively; see Fig.|7]i. In ad- 
dition, spectra were extracted in three other important epochs: 
a) around the peak of the 7-ray emission (see Fig.|8]l, yielding a 
spectrum with a threshold of 400 GeV (T4 0-Peak); b) in the 
central five hours characterized by a threshold as low as 200 GeV 
(T2 dataset); c) in the epoch simultaneous with the RXTE ex- 
posure (T300-RXTE window, 44 minutes overall, see Table 1), 
where the combined X-ray spectrum can be measured over 2 
decades in energy. 

The results of the spectral fits are given in Table |2] with a 
selection shown in Fig.|9] The spectra present a significant cur- 
vature with respect to the pure power-law. The probabilities 
show that the latter is completely excluded in the high states, and 
is unlikely in the low state. The spectral curvature is generally 
well described either by a power-law model with an exponen- 
tial cutoff around 1 TeV ( ^(E) = $0 ^'"^ e--^/-^<^"' ), or a 
log-paraboHc function ( <^{E) = $0 E'^'^ + ^ ). Most 

remarkably, the curvature of the spectrum is strongly variable 
with time. In particular, the curvature is more pronounced (i.e.. 



the parameter h is larger) in the brightest state and decreases as 
the source dims. This represents direct proof that the curvature 
of the VHE spectrum in PKS 2155-304 is also of intrinsic ori- 
gin, inside the emitting region, and cannot be attributed entirely 
to 7-7 absorption on the EBL or on any local external field that 
is constant on the observed timescales. 



Besides providing the widest energy coverage, the T2 
spectrum allows a direct comparison with the spectrum mea- 
sured during the first exceptional flare on the night of July 27- 
28. The latter has the same energy threshold (200 GeV) and is 
well described by a broken power-law ( $(£^) = $0 E^^^ for 

E < Sbrcak and = $0 E^^'I'E-^^ for E > Streak ) 

withTi = 2.71±0.06,r2 = 3.53±0.05, and £:brcak = 430±22 
GeV. Fitting this function to the T2 00 spectrum yields almost 
identical results, of the same break energy, slopes and change 
in spectral index by AF ~ 0.9 (Table |2]i. This shows that the 
source was in a similar state, even though the overall average 
normalization is about ^30% lower than two nights before. As 
for the July 27-28 night, the T2 00 spectrum is not well fitted 
by a power-law model with exponential cutoff or a log-parabolic 
function (F-test >99% compared to the broken power-law). Both 
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Table 3. Spectral fit of the same VHE spectra given in Table|2l but corrected for EBL absorption (statistical errors only). 
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functions underestimate significantly the 7-ray flux at higher en- 
ergies. 

4.1.1. Correction for intergalactic 7-7 absorption 

The VHE 7-ray emission from extragalactic sources is expected 
to be attenuated by photon-photon interactions with the EBL 
photons in the optical-to-IR waveband. The energy dependence 
of the optical depth - which is determined by the spectrum of the 
diffuse background - causes a general steepening of the emitted 
7-ray spectrum, more or less severe acc ording to the energy band 
considered (see e.g.. lAharonianllToOlL and references therein). 
At redshift z=0.116, this effect is substantial and must be taken 
into account to study the true energy output, sp ectral properties, 
and location o f the Compton peak in the SED dAharonian et alj 
I2005bll2006ah . Source diagnostic based on flux or spectral vari- 
ability, instead, is unaffected, since intergalactic absorption is a 
constant factor for all purposes related to blazar variability (the 
diffuse background varies only on cosmological timescales). 

The EBL waveband that affects the observed VHE band the 
most is dominated by the direct starlight emission. To correct for 
EBL absorption, as reference we adopted the model of t he EBL 
spectral energy distribution bv lFranceschini et"aLl (|2008|) . which 
is based on the emission from galaxies. This model takes into 
account the most recent results on galaxy properties and evo- 
lution and is consistent with bo th the low er limits from source 
counts - in the UV-optical ( Madau & Pozzettj 2000,) as well as 



near-mid infrared waveband dPazio et al .120041-lDole et al.l2006l) 
- and with the up per limits derived from the TeV sp ectra of high- 
redshift blazars ( (Aharonian e t alj |2006al l2007bllH l2()02al) It is 
similar in shape to both the model by 'Primac k et al. ( 20051) and 
the "low-IR" calculation by .Kneiske et al.. (.20041) . The spectra 
were corrected by applying the optical depth calculated for the 
average observed photon energy in each energy bin. 

However, it is important to recall that a significant uncer- 
tainty in the SED of the EBL still remains, since it could be 
both lower and higher than assumed: either down to the abso- 
lute lo wer limits given by HST galaxy counts ( Madau & Pozzettil 
I2OOOI) (as in the model bylMmack et al."2005^ or up to the up- 
per limits given by TeV blazars (tAharonian et al, 2006a). To es- 
timate this uncertainty in both s hape and normaUzatio n, we also 
used the shape of the model by iPrimack et al.l (l2005h . rescaled 



to these two levels. In the energy range around the starlight 
peak (1-3 fim ), the residual uncertainty in the EBL absolute 
normalization is of the order of 50% (from ~8 to ~12 nW 
m~^sr~^at 2.2 micron, while our reference model gives 9.4 nW 
m~^sr~^). This translates into a systematic uncertainty of the or- 
der of Ar ~ ±0.2 in the reconstructed 7-ray spectrum. Namely, 
the reconstructed spectra (which we call "intrinsic") discussed 
in the following Sects, can actually be up to ~0.2 steeper or 
harder than indicated. When relevant, we take this systematic 
uncertainty into consideration, but as we show in the following, 
it does not change the main properties of the 7-ray spectrum and 
Compton peak frequency of PKS 2155-304. 

4.1 .2. Absorption-corrected 7-ray spectra 

For the power-law spectra measured on short timescales (as 
given in Fig. |2l), the intrinsic spectra are again well fitted by a 
power-law model with a slope that is typically harder by —0.8 
(namely. Tint = Fobs — 0.8). The results of the fits to the spectra 
with higher event statistics are provided in Table [3] Even after 
correction for the steepening induced by EBL absorption, the 
7-ray spectra show clear evidence of curvature, and the power- 
law model is excluded with high confidence for all states but the 
RXTE epoch (which has the lowest exposure). The spectral cur- 
vature is well described by either the log-parabolic function or a 
power-law model with an exponential cutoff around 1-2 TeV, for 
all spectra, including the T2 00 dataset. 

When a spectrum shows a relatively uniform curvature as in 
this case, however, the log-parabolic model is generally prefer- 
able. It has the advantage of providing a more direct measure of 
the curvature in the true observed band, whereas the exponential 
cutoff model tends to match a given curvature in the observed 
passband by using a specific section of its cutoff region, and 
pushing the power-law component outside the actual observed 
range. This often yields artificial values for the slope, which are 
typically too hard. The log-parabolic fit allows also a straightfor- 
ward estimate of the location of the SED peak (Epeak, defined by 
r(Epeafe) = 2) from the curvature itself, with a minimum of free 
parameters. To this aim we used the functional form described in 
Tramacere et al. (2007), where b and Eppak are the independent 
free parameters instead of b and FitoV- The comparison of the 
curvatures among different states and between the synchrotron 
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Fig. 10. Evolution of the X-ray spectrum with time, for the whole Chandra exposure. Left panel: each 8-min bin is fitted with a 
single power-law model plus galactic absorption. The HESS window ends at MJD*=46. 16. Left lower panel: the colors and symbols 
(in time sequence: red circles, magenta triangles, green squares, cyan triangles, blue circles) mark the different zones for the index- 
flux correlation shown in the right panel. Right panel: for clearer visibility, the last two intervals (cyan triangles and blue circles) 
are further binned into 16-minute spectra. There is a clear "harder when brighter" trend in the decaying phase of the first branch 
(corresponding to the VHE window). The behaviour changes in the successive rising phase: the X-ray spectrum continues to soften 
while the flux increases, drawing a counter-clockwise pattern. 



and IC components can also provide important clues about the 
source emission regime (Thomson or KN) and the acceleration 
mechanisms (Massaro et al. 2006). 

By comparing the spectra in the 3 different flux states 
(T400-Peak, T300-High and T300-Low), one can see that 
the curvature changes significantly along the night (at a confi- 
dence level > 99.99%) and so does the Compton peak energy. A 
clear trend emerges: both the curvature parameter h and Epcak in- 
crease with the VHE flux. At the maximum of the VHE flare, the 
spectrum is strongly curved (5 = 1.2 ± 0.2), with the Compton 
peak estimate at Epeafc = 500 ± 50 GeV. As the flux decreases, 
the curvature flattens (h =0.62 to 0.35), while the IC peak shifts 
to lower energies (Epeafe = 260 ± 35 GeV to Epeafe = 70 ± 50 
GeV, respectively). A lower/higher EBL level does not change 
these results substantially: a higher level yields similar curva- 
tures (6=1.4, 0.76, and 0.48, respectively) and slightly higher IC 
peak energies because of the generally harder spectra (Epeafc= 
580, 370, and 180 GeV, respectively). It is important to recall 
that the absolute value of the curvature b depends on the par- 
ticular choice of the EBL spectrum used, but not the trend itself. 
This trend is opposite to what is generally observed and expected 
for the synchrotron emission in TeV blazars. For example, in 
Mkn421 the curvature h dec reases as both Epeafc and the flux 
increase ( iMassaro et alj|2064h . We also note that the spectral in- 
dex r does not correlate (and possibly anti-correlate) with the 
curvature h, in contrast to what is ob served in the X-ray b and for 
this (see next Sect.) and other HBL (iMassaro et alj|2004l) . 

4.2. X-ray spectra 

The Chandra spectra were extracted both by a uniform sampling 
of the whole exposure, in different time bins (2-4-8-16 minutes). 



and strictly coincident with the VHE time bins. The spectra were 
all fitted using an equivalent hydrogen column density fixed at 
Galactic values (1.69 x 10^"cm^^) with different source models. 
On short integration times (< 1 hr), a single power-law model 
provides statistically good fits for all datasets, while evidence of 
curvature is found only when larger exposures are considered 
or a wider energy band is available (for example including the 
RXTE data). 

The time evolution in the X-ray spectrum during the entire 
Chandra pointing -which extends few hours beyond the end of 
the HESS observation- is shown in Fig. [10] There is a clear 
trend of "harder-when-brighter" behaviour in the first part of the 
dataset, corresponding to the decaying phase of the main VHE 
flare. This behaviour is also followed by the small-amplitude 
flares, whose paths in the flux-index plane overlap tightly with 
the overall trend of the decaying phase (Fig. [TOl right panel). 

However, the relation changes in the last part of the ob- 
servation (MJD*>46.16): as the X-ray flux starts to increase 
again, the spectral index continues to soften. This "softer-when- 
brighter" behaviour in the rising phase of a new flare reveals a 
change in conditions for the emitting region. It is indicative of 
a slow acceleration/injection process, whose timescale is com- 
parable with the other timescales of the system {tacc ~ tcooi)- 
The information about the flare then propagates fro m lower to 
higher energies as particles are gradually accelerated dKirk et alj 
119981: iRavasio et al]|2004h . If the optical variations are indeed 
associated with the flaring zone, the optical data would support 
this scenario as well, exhibiting increasing flux just before the 
X-ray rise at the end of the Chandra observation. Together, the 
two patterns of the X-ray data draw part of a counter-clockwise 
loop in the flux-index plane (Fig.fTOli. 
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Fig. 11. Left panel: Chandra spectrum for the total exposure, fitted with a broken power-law plus galactic column density. The 
additional absorption feature in the 0.3-0.4 keV range is due to the contaminants on the ACIS optical blocking filter, not yet fully 
accounted by the calibration. It can be accounted for with a simple edge model at 0.3 1 keV and rniax=0.4 (see text). Right panel: the 
simultaneous Chandra+RXTE spectrum in the common 44-min window. The plot corresponds to a single power-law with galactic 
absorption, and the data/model ratio shows the clear evidence of curvature (fit parameters are given in TableH]). The RXTE /Chandra 
normalization is fixed at 1.08, as derived from the fit in the overlapping energy range (3-7 keV). 



Table 4. Fit of the X-ray spectra simultaneous to the VHE data. 
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Fig.[TT]shows the spectra of both the total Chandra exposure 
and the RXTE simultaneous window. The results of the fits per- 
formed on the HESS -simultaneous datasets are given in TablelH 
All datasets correspond to strictly simultaneous windows except 
for the T300-High spectra, for which the X-ray window does 
not include the first ~10 minutes of the 1.3-hrs VHE window 
above 300 GeV. Since there are no significant spectral changes 
at VHE in that window, the VHE spectrum can be considered to 
accurately represent the 7-ray spectral shape in the X-ray win- 
dow. 

For all spectra in TablelH there is clear evidence of curvature, 
and the single power-law model is rejected with high confidence 
(F-test> 99.99%). The spectra show a continuous steepening to- 
wards higher energies up to ^^20 keV, which is well represented 
by both a broken power-law and log-parabolic models. A power- 
law with an exponential cutoff is excluded as well (P^2 <0.009) 
for the spectra with the highest statistics (T2 00 and total spec- 
trum). The drop rate in the cutoff region is significantly slower 
than Q-^/^^'t and also slightly slower than Qri^/^<:^t) ' , as 
indeed expected for the synchrotro n emission of a p article dis- 
tribution with an exponential cutoff d Aharonianll2000l) . 

As is clear from Fig. [TO] and Table |4l during all times the 
X-ray spectrum of PKS 2155-304 remains steep, with a con- 



vex shape and no signs of flattening at high energies (as instead 
found in XMM-Newton observations pe rformed in November 
2006: iFoschini et allllOOl: IZhang||2008h . This means that the 
peak of the synchrotron emission has not entered the observed 
energy range at any time, and that there is no sign of the possible 
emergence of the IC component in the hard X-ray band. 

It is interesting to compare the curvature parameters and 
SED peak location given by the log-parabolic fits. Both the 
spectral index and the curvature increases (slightly) as the flux 
decreases. This is also corroborated by the fits of spectra ex- 
tracted in even shorter intervals at the two extreme of the X-ray 
flux range (namely T300-Xmax and T400-Xmin). The log- 
paraboUc fit yields F = 2.45 ± 0.02 and 6 = 0.16 ± 0.04 ver- 
sus F 2.64 ± 0.01 and h = 0.24 ± 0.04, respectively, for 
an integrated flux Fo.5_5koV =4.33 and 1.65x10"^" erg cm^^ 
s~^. The change in the two spectral parameters, however, is such 
that the estimate of the location of the SED peak remains basi- 
cally constant: for all spectra, the synchrotron Epeak falls within 
the range 40-50 eV (with a typical l-cr statistical error of ±20 
eV). In contrast to the behaviour in 7-rays, in the X-ray band 
the photon index shows a positive correlation wit h the curvature 
6, as typically observed for example in Mkn 421 dMassaro et al.l 
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Fig. 12. Plot of the simultaneous spectral and flux variability in 
the X-ray (blue circles) and VHE (red triangles) bands, in the 
T300-X window. The spectra have been extracted in the same 
time -bins, with integration times of 7 and 14 minutes (before and 
after MJD*46.0, respectively). The vertical scales are different 
(left: X-ray values; right: VHE values). Upper panel: the scale of 
the VHE photon indices is 2.6 times the X-ray scale. The VHE 
indices refer to the observed spectra; the corresponding values 
after correction for EBL absorption can be obtained as Tintr = 
TyHE — 0.8. Lower panel: integrated energy fluxes in the same 
units of erg cm~^ s~^. The range of the VHE scale is the cube 
of the X-ray range (15.3 x versus 2.5 x). The vertical lines mark 
the three time-zones (A, B, and C) referred to in the text and in 

Fig.ini 



|2004|) . The absolute values of the curvature are similar to those 
found for most other HBL (iMassaro et al.ii2008i) . 

5. X-ray vs. TeV correlations 

5.1. Spectral variability 

The exceptional 7-ray brightness observed in this night, coupled 
with the sensitivity of the HESS array and the continuous cover- 
age provided by Cliandra, allows the emission in the two bands 
to be compared with unprecedented time resolution in the spec- 
tral domain as well. 

Fig- m shows the simultaneous flux and spectral properties 
measured in 7 and 14 minute time bins, in the T300-X win- 
dow. Inside this window, the VHE spectral index can be well 
constrained (±0.1) over approximately a decade in energy (0.3- 
2 TeV). The different binning was chosen to achieve compara- 
ble S/N ratio during the night, and as a good compromise be- 
tween spectral determination at VHE and time resolution (see 
Sect. 4.1). Both X-ray and VHE spectra have been extracted in 
exactly the same time bins. The spectra were fitted with a single 
power-law model, which provides a good fit for each time bin. 
The integrated energy fluxes were calculated using the specific 
spectral value measured in each bin. 

Two general properties can immediately be noted. The first 
is that the VHE emission shows a definite correlation with the 
X-ray emission not only in flux but also spectrally. The spec- 
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Fig. 13. Plot of the 7-ray vs. X-ray flux correlation, in a log-log 
diagram and units of erg cm^^ s^^. Different markers corre- 
spond to the three different intervals shown in Fig. [12] A (black 
circles), B (blue triangles) and C (green square). For visual refer- 
ence, the dotted lines show three different slopes of the relation 
F~f cx F^. For clearer visibility, the Y-scale corresponds to the 
square of the X-scale, so that a quadratic relation has the slope 
of 1. The best-fit values are given in Table |5] 



tral evolution follows the same overall pattern in the two bands, 
although with different amplitudes. The correlation coefficient 
between the X-ray and 7-ray spectra is r = 0.65, with a proba- 
bility P < 0.1% of a chance correlation. 

The second property, as previously illustrated by Fig. [3] is 
that the source shows amplitude variations much larger in 7-rays 
than in X-rays. This is now evident also for the spectra, although 
not as dramatically as for the flux: the spectral variation at VHE 
is about 3 times the variations in the X-ray band ( AFvhe ~ 0.65 
vs. AFx w 0.21). 

The combination of such large-amplitude variability at VHE 
with correlated but small-amplitude variations in X-rays yields 
one of the most striking features of this dataset: the VHE flux 
varies more tlian quadratically with respect to the X-ray flux 
during a decaying pliase. 

5.2. Cubic relation between X-ray and TeV fluxes 

Fig. [T3] shows the VHE flux as a function of the X-ray flux, 
in a log-log diagram. The data are divided into three subsets 
("A", "B", and "C") corresponding to three characteristic epochs: 
epoch 'A" covers the first 35 minutes of the simultaneous win- 
dow; epoch "B" corresponds to the full rise and decay phases of 
the small-amplitude flare at MJD*45.96; epoch "C" covers the 
remainder of the dataset, corresponding to the broader flare (see 
Fig. [Ell. 

The data were fitted with a linear relationship in the log-log 
space (F^ cx F-^), for the total dataset and in each subset sepa- 
rately. The results are reported in Table |5] together with the fits 
of the same datasets with a finer sampling (4-minute bins). 

The 7-ray flux traces the variations in the X-ray flux far more 
than quadratically, namely approximately as Fvhe F^. With 
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Table 5. Values of the slope j3 of the correlation cx . 
Fluxes integrated over each respective energy band (0.5-5 keV 
and 0.3-3 TeV) and in strictly simultaneous bins. The parameter 
is shown for fits performed with two different binnings (4 min 
and 7-14 min) and over three different intervals, as shown in Fia. 

m 



Datasets 


4-min bins 


7-14min bins 


all 


2.21 ±0.05 


2.25 ±0.05 


A 


no corr. 


no corr 


B 


2.72 ±0.17 


3.18 ±0.18 


C 


2.83 ±0.17 


3.14 ±0.18 


B+C 


3.13 ±0.11 


3.35 ±0.11 



the exception of the first few minutes of epoch A, all points lie on 
a narrow path, also during the small flare on MJD*45.96 (epoch 
B). In epoch A, instead, the two emissions do not seem to cor- 
relate. The lack of correlation is limited to a few points and in a 
short interval; however, one may speculate whether this is caused 
by the rapid, low-significance structures in the VHE light curve, 
which are not evident in X-rays (as discussed in Sect. 3.3). 

One of these subflares is indeed present in epoch A around 
MJDo = 45.925 (see Fig. |6]l, and its exclusion does bring the 
7-ray emission more in line with the X-ray pattern. This type 
of 7-ray subflares o n top of correlat e d emi ssions have been re- 
cently envisag ed by iGhisenini etal] (|2009|) . but more sensitive 
instruments are required to draw any conclusions. 

In the epochs B and C, the cubic correlation is obtained 
by considering both the two zones separately and together, for 
which an even steeper slope of /? = 3.35 is obtained. This cubic 
correlation is robust with respect to the inclusion or exclusion 
of single data points and, in particular, does not depend on the 
lowest VHE point. 

The flux-flux correlation is plotted using the same units for 
both X-ray and VHE bands, namely the integrated energy fluxes. 
Compared to previous studies, which used the observed event 
rates, our approach is more consistent and it is allowed by the 
good spectral determination in each time bin. However, we re- 
mark that the measured relation does not depend significantly 
on th e particular approach used (see discussion in iFossati et alj 
I2008h : a cubic relation is also obtained by using the event rates, 
both in the VHE band (as photons cm~^ s~^) and in the X-ray 
band (count rate). 

The correlation found in the total dataset is less steep than 
in each subset, while in the Bh-C epoch it is slightly steeper than 
for each of B and C separately. This can be caused by a possi- 
ble shift in the flux-flux paths among different datasets/epochs. 
While individual paths still obey a specific steep trend, taken to- 
gether they can produce a flatter (or s teeper) envelope. T his ef- 
fect was indeed observed in Mkn 421 ( iFossati et alj|2008l) . when 
considering data for different days. 

During the decaying phase of a flare, a cubic relation be- 
tween the 7-ray and X-ray fluxes is not easy to explain even for a 
source in Thomson condition, if b oth fluxes sample the emi ssion 
beyond the respective SED peaks jKatarzvnski et alj|2005h . It is 
the first time that such a steep slope has been observed in the 
history of the X-ray/TeV correlation studies, though indication 
of a super-quadratic relation was recently reported for Mkn 42 1 
( IFossati et alj|2008l) . during single flares. 

These correlations have been studied so far mainly with the 
RXTE -PCA instrument, which samples higher energies than 
those observed here. If the spectrum changes with the flux, and 
with a "harder-when-brighter" behaviour as in this case, the am- 



Table 6. Values of the slope (3 of the correlation F^ cx F^, 
for fluxes integrated over different energy bands, as indicated. 
The 3-15 keV flux is obtained by extrapolation of the power-law 
model in the Chandra passband. Fit of the dataset "Bh-C" in the 
7-14 minutes binning. 



VHE bands 


X-ray bands 


0.3-0.7 TeV >0.7 TeV 


0.5-5 keV 
3-15 keV 


2.91 ±0.12 4.11 ±0.27 
1.87 ±0.08 2.70 ±0.20 



plitude of the variations changes with energy, and thus the slope 
of the correlation can depend on the observed band. 

To quantify this effect, we also investigated the flux-flux re- 
lation by extrapolating the Chandra spectrum in the RXTE band 
(namely, integrating the best-fit model in the 3-15 keV range), 
and dividing the VHE range into soft and hard bands (0.3-0.7 
and >0.7 TeV, respectively). The result is shown in Table|6l The 
same cubic correlation observed in Chandra and HESS trans- 
lates into a quadratic relation (in fact quite similar to that ob- 
served in Mkn 42 1 ) between the hard X-ray and soft VHE bands, 
while an even steeper slope {(3 > 4) is obtained between the soft 
X-ray band and hard VHE band. We note however that the cubic 
relation is not simply the effect of a pivoting type of variabil- 
ity seen in different energy bands. The VHE band is closer to 
the IC peak (i.e the pivoting point) than the X-ray band to the 
synchrotron one, thus the X-ray variations should be larger than 
the VHE ones, in contrast to what is observed. It is the entire 
IC peak that has actually varied far more than the synchrotron 
peak, as shown also by the nearly constant and cubic values of 
the correlation between corresponding bands (i.e., soft-soft and 
hard-hard, see Table 3). 

Comparing thes e 2006 data ( low state) with the X-ray/TeV 
campaign of 2003 dAharonian et al. 2005b), the overall bright- 
ening in the X-ray (2-10 keV) and VHE (>300 GeV) bands is 
similar (a factor of 2.7 and 5, respectively), and corresponds to 
a relation F^ cx F^-^ between the two epochs. Thus PKS 2155- 
304 has varied its overall synchrotron and IC luminosity sub- 
quadratically on very long timescales, but super-quadratically on 
intra-night timescales, at least during this major flaring event. 



6. Spectral Energy Distributions 

6.1. X-ray/j- ray spectra pairs 

To highlight the evolution of the SED during the night, a gallery 
of selected pairings of simultaneous X-ray and VHE spectra is 
shown in Fig. [14] The scales on both axes are kept the same for 
both X-ray and 7-ray energies, to enable visually a correct com- 
parison of the spectral slopes. For reference, a selection of his- 
torical observations is also plotted, in particular the data from 
the first X-ray/TeV multiwavel ength campaign perfor med on 
PKS 2155-203 in October 2003 dAharonian et al.ll2005bh . These 
data (obtained with RXTE and HESS) correspond to one of the 
historically lowest states ever observed from this object. 

The panels show the SED-snapshots taken in the bright- 
est/hardest state (T300-Xmax), in the faintest/softest state 
(T400-Xmin), and simultaneously with the RXTE pointing 
(T300-RXTE). In addition, Fig . [T4]also shows theT400-Peak 
spectrum (i.e., extracted around the maximum of the 7-ray flare; 
see Tableland Fig.|8]i. This spectrum unfortunately lacks X-ray 
coverage (the Chandra pointing started a few minutes later), but 
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Fig. 14. Gallery of simultaneous X-ray and VHE SED-pairs in different states along the night. Left panel; maximum (T 3 -Xmax 
blue triangles) and minimum (T4 00-Xmin, red squares) states (see Table 1 and Fig. 2). Right panel: simultaneous Chandra+RXTE 
and HESS spectra (T300-RXTE, blue squares). For comparison, the following data are also plotted: the VHE spectrum at the peak 
of the 7-ray flare (T4 0-Peak; upper black squares); the multiwavelength R XTE -HESS campaign i n 2003 (lower black circles; 
lAharonian et al.ll2005bl) ; historical EUVE flu x and BeppoSAX spectrum (from lChiappetti et alj[l999l grey points). AU VHE data 
are corrected for EBL absorption according to lFranceschini et alT(l2008l) (see Sect. 4.1.2). 



it corresponds to the brightest 7-ray emission ever recorded from 
PKS 2155-304. 

Both synchrotron and IC peaks do not shift across the ob- 
served passbands, despite the large flux variations. There is 
no evidence of the dramatic changes displayed by Mkn501 or 
lES 1959H-650. Only at the flare maximum does the Compton 
peak become visible in the observed passband (between ^400 
and ^^600 GeV, depending on the EBL normalization). The X- 
ray spectrum hardens apparently pivoting around the UV band. 
This behaviour is typically observed also in Mkn421, Mkn501 
and lES 1959H-650, but in these sources the amplitude of the 
spectral hardening tends to be significantly more pronounced, 
leading to a shift of the synchrotron peak in the hard X-ray band. 
In the VHE band instead the behaviour is more complex, since 
the hardening at higher fluxes is also accompanied by a stronger 
curvature/cutoff (see Fig. [141 right panel). 

The slopes of the X-ray (above ~1 keV) and VHE spec- 
tra are very similar, going from F ~2.6 (± < 0.04) to 2.9 
(± < 0.11) between the high and low states. It is also interesting 
to note that the X-ray spectrum is significantly less curved than 
the sp ectrum measured by B eppoSAX during the high state of 
1998 (IChiappetti et al.lll999l) . While the slopes above few keV 
are quite similar (see Fig.O, the flare in this night seems char- 
acterized by a higher luminosity in the soft X-ray band. 

6.2. Compton dominance 

The time evolution of the lyF^ fluxes close to the SED peaks is 
provided by Fig. [3] As one can see, the Lc/Lg ratio is of the 
order of ^8, but it is rapidly variable - on the same timescales 
as the flux variations - decreasing to the usual values of < 1 in a 
few hours. 

We estimated the Compton dominance using also the more 
detailed spectral shape obtained for the T300-High dataset, 
and renormalizing the flux to the 2 and 4-minute light curves. 
Considering the integrated luminosity over a decade in energy, 
namely 0.3-3 TeV and 0.3-3 keV, the ivHs/ix ratio is wlO 
in the first minutes of the simultaneous window. Assuming that 
the cubic trend between 7-ray and X-ray fluxes is maintained up 



to the flare maximum (which occurred ^^25 minutes before the 
start of the X-ray observation), one obtains LvHE/ix~14 at the 
flare maximum. The lyF^ flux at the peak of the SED can instead 
be estimated from the log-parabolic fit, based on the assumption 
that the curvature of the spectrum also remains the same outside 
the observed energy band. Using the aforementioned procedure, 
one obtains a LqILo, ~5 in the simultaneous window, increasing 
to ^^6 at the flare maximum. The values are lower because the 
estimate of the synchrotron peak energy (^45 eV) locates the X- 
ray band farther away from the synchrotron peak than the VHE 
band is from the Compton peak (--400-600 GeV). 

This is the first time that such high LqIL^ ratios are observed 
in an HBL, irrespective of the amount of the intergalactic EBL 
absorption. While the Compton dominance can be up to 100 in 
powerful FSRQ, boosted by the IC emission o n the intense exter- 
nal photon fields from the disk and the BLR jSikora et al.llI994l; 
iGhisellini et al.in"998h . so far it has been of the order of unity 
or less in HBL. High Compton dominances in HBL were pre- 
viously o btained only in the presence of a very high density of 
the EBL ( lAharonian et al.ll2002al l2006ah. which is now consid- 
ere d very unlikely (lAharonian et al .ll2006aHMadau & Silkll2005l: 



iFranceschini et al.ll2008h 



6.3. Overall SED properties 

The overafl SED of PKS 2155-304 in the highest and lowest 
states during this night is shown in Fig. [15] together with histor- 
ical data. To plot the 7-ray spectral shape in more detail while 
preserving the amplitude variation, the VHE spectra shown in 
Fig.[T5]are the T30 0-High and T30 0-Low average spectra de- 
scribed in Sect 4.1 (see Fig.|2]and|9]l, rescaled to match the high- 
est and lowest fluxes in the 4-minute light curve. Since within 
these two subsets, the spectra are compatible with a constant 
value (see Fig. |7]i, this procedure should not introduce signifi- 
cant distortions in the VHE spectral shape. For clearer visibil- 
ity, the SED focuses on data above ^10^^ Hz, thus excluding 
the radio frequencies (VLBI range). With the flux and variabil- 
ity timescale (<3 hrs) shown at optical (and higher) frequencies, 
the synchrotron emission coming from this region becomes self- 
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absorbed already at frequencies below ^^10 Hz, and is thus not 
expected to contribute significantly in the radio range. 

From the results of the log-parabolic fits, it is possible to esti- 
mate the peak frequency and luminosity of both the synchrotron 
and Compton emission. By considering the highest state, the 
synchrotron and Compton peak luminosities are estimated to be 
at 1.5 X lO**^ and 7x 10**^ erg/s, respectively, with peak frequen- 
cies of ^50 eV and ~260 GeV. Extrapolating the X-ray data 
with the observed trends up to the flare maximum, one obtains 
1.6x lO-^s and 9x10^^ erg/s at ~50 eV and -500 GeV, respec- 
tively. In the lowest state, instead, the peak luminosities become 
comparable (1.2x lO"*^ and 1.1 x 10*^ erg/s, respectively, at —50 
eV and —70 GeV). As discussed in Sect. 4.1, even the lowest 
possible EBL density does not alter substantially these estimates. 

In the optical band, the contribution of the host galaxy is neg- 
ligible; thehost galaxy i s resolved in optical (F alomo et al. 199 1]; 
lFalomolll996 ') and NIR jKotilainen et al.ll 1998b , and found to be 
a giant elliptical galaxy of M(R)=-24.4. This translates into an 
apparent m(V)~ 15.7 — 15.8 u sing the typical colou rs for an 
efliptical galaxy ( V-R=0. 6 1-0.71, 'Fukugi ta et al.ll 1 9951) . The op- 
tical flux of PKS 2155-304 is thus always dominated by the jet 
emission (see e.g.. lDolcini et al.l2007l) : the historical light curves 
from long-term photometric monitoring i n the V band show vari- 
ations in the range 12 .3-13.9 magnitudes dCarini & Milled 1992t 
lOsterman et al.ll2007l) . During this night, the optical flux is very 
high but still far from the highest fluxes observed from this ob- 
ject. 

In the X-ray band, instead, the flux is close to the high- 
est state observed historically, as measured by RXTE -PCA i n 
1996 {F2-1Q -1.6x10-1" erg cm-^ s^^ lUrrv et all Il999l) . 
PKS 2155-304 is observed frequently in the X-ray band, since 
it is a calibration and monitoring source for XMM-Newton, 
Chandra, and Swift, but was never found at the level observed 
durin g this night, in the 2- 10 keV band (see e.g. Donato et aL 
I2OOII: iMassaro et"!!] |2008|) . Quite interestingly, the Chandra 
spectrum in the highest state also seems to connect smoothly 
with the flux and spectrum measured by RXTE -HEXTE in 1996 
dUrrv et al.ll 19991 butterfly in Fig. [TSl l. It is also interesting to 
note that the extrapolation to lower energies of the log-parabolic 
fit for the average T2 spectrum matches quite well the flux in 
the optical band. 



7. Summary of the main observational findings 

Before discussing the implications in the context of blazar 
physics, it is useful to summarize the main observational find- 
ings of this phenomenologically rich dataset. 

- Large-amplitude 7-ray variations are accompanied by small 
X-ray and optical changes. In a few hours, the 7-ray flux 
changes by more than an order of magnitude, reaching a lu- 
minosity of « 10*^ erg/s, while the X-ray flux varies by only 
a factor —2 overall, and the optical V flux by less than 15%. 

- The X-ray and 7-ray emission correlate strongly, overall and 
on short (sub-hour) timescales. On very short timescales 
(few minutes), the behaviour might be more complex. 

- There is no evidence of time lags between the X-ray and 7- 
ray emissions, with a 95% upper limit of —3 minutes for the 
overall light curve. In addition, no lags are found between 
the hard and soft energies of each passband. 

- The optical light curve shows a —15% rise that appears to 
start simultaneously with the 7-ray flare, but develops on 
much longer timescales, reaching a plateau about 2 hours 



later than the VHE peak. The optical emission does not show 
any correlation with the other two bands on short timescales. 

- When correlated, the "f-rayflux decreases as the cube of the 
X-ray flux (Fvhe « F^^). This cubic relation holds both 
during the overall decaying phase and considering shorter 
intervals separately. 

- The X-ray and 7-ray (> 300 GeV) spectra correlate as well, 
following similar patterns in their time evolution, but again 
of different amplitudes. The VHE variations (AFvhe = 
0.65) are wider by a factor of 3 than those in the X-ray band 
(AFx = 0.21). 

- The 7-ray spectra are significantly curved, and the curva- 
ture changes with time, in correlation with the flux state: the 
higher the flux, the more curved the spectrum (the curva- 
ture parameter h goes from 1.2 ± 0.2 at the maximum to 
0.38 ± 0.14 in the low state). 

- The synchrotron and Compton peaks show no strong shift 
in frequency across the observed bands, despite the dramatic 
luminosity changes (lOOx the quiescent state, 20x in this 
single night), remaining close to their historical values. From 
the curvature of the spectra, the synchrotron peak can be es- 
timated at « 40 — 50 eV constantly along the night, while 
the Compton peak shifts from w 500 ±100 GeV at the max- 
imum to « 70^10° at the end of the night (uncertainties 
given by the range on the EBL normalization). 

- A very large Compton dominance is observed (Lc/^sJ^8). 
This is the first time that such a high Lq/Lo, ratio is seen 
in an HBL, irrespective of the level of intergalactic EBL ab- 
sorption. However, it also evolves rapidly, decreasing in a 
few hours to the more usual values of < 1 . 

- X-ray and VHE spectra shows a similar "harder-when- 
brighter" behaviour in the simultaneous window. The X-ray 
data alone also sample the start of another flare, character- 
ized instead by a "softer when brighter" behaviour 

8. Discussion 

The results obtained from this campaign seem to both corrobo- 
rate and challenge the one-zone SSC interpretation at the same 
time. 

On the one hand, the strong correlation between variations in 
the X-ray and VHE bands - the emission follow the same vari- 
ability patterns in terms of both flux and spectrum, and without 
apparent lags - do indicate that the same particle distribution, in 
the same physical region, is likely responsible for the activity in 
both energy bands. The simultaneity of the occurrence of the op- 
tical and VHE flares also suggests that the emission in all three 
bands is responding to the same flaring event. 

On the other hand, the cubic relation during the decaying 
phase cannot be easily accounted for within a one-zone SSC sce - 
nario. This is shown for example by iKatarzvriski et alj ( l2005h . 
who investigated the different X-ray/VHE correlations achiev- 
able in a one-zone SSC scenario, for a plausible range of phys- 
ical parameters and considering all possible combinations of 
fluxes from both sides of the SED peaks. The relation between 
X-ray and 7-ray fluxes when both bands sample the emission 
above their respective SED peaks is typically less than quadratic. 
Even imposing the Thomson condition for the IC scattering of 
X-ray photons, which in this case would require the extreme 
values of 5 —800 and B —0.6 mG (see next par.), at most a 
quadratic relationship can be obtained, but not in the decaying 
phase if the latter is due to radiative cooling. In general, corre- 
lations steeper than quadratic (cubic or even more) are indeed 
possible, but only when the X-ray band corresponds to frequen- 
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Fig. 15. Synthetic SED of PKS 2155-304 showing the high- 
est and lowest simultaneous states duri ng this night, together 
with historica l data ( shown in grey, se e Chiappetti et al.|[l999l : 
iBertone et alj I2OOOI : lAharonian et alj I2005bl and references 
therein). The hard X-ray data (butterfly) correspo nds to the 
RXTE -HEXTE spectrum in the high state of 1996 dUrrv et al.l 
119991) . The right axis gives the luminosity scale in erg/s. Highest 
state (blue triangles): T300-High spectrum scaled to the high- 
est VHE flux in the 4-min light curve. Lowest state (red squares): 
T3 0-Low spectrum scaled to the lowest VHE flux in the 4-min 
light curve. The X-ray states in the corresponding time bins are 
practically equal to the T300-Xmax and T400-Xmin spectra, 
which are thus plotted here. The dashed line shows the one-zone 
SSC fit of the 2003 m ultiwavelength campaign (black circles, 
lAharonian e"td]|2005bh . 



cies below or close to the synch rotron peak, as in the c ase of the 
1997 flare of Mkn 501 (see e.g.. Pfavecchio et al.ir200lh . 

A further ingredient is therefore needed to obtain a cubic 
relation, in a one-zone SSC scenario. It is important to remind 
here that "one-zone" does not mean the request that the entire 
SED is produced by a single zone. It is generally understood 
that the SED is always formed by the superposition of radiation 
from different zones along the jet, both in space (e.g., the radio 
emission has to come from a much larger region than for the 
rapidly variable X-ray or optical emission, to avoid synchrotron 
self-absorption) and in time (multiple injections). What is gen- 
erally tested with a "one-zone scenario" is the hypothesis that at 
any given time - but during large flares in particular - it is one 
single region or population that dominates the entire radiative 
output, determining both peaks of the SED. This scenario has 
been very successful in explaining so far the spectral and tempo- 
ral properties of blazars, both among different objects and in sin- 
gle sources during large flares, (e.g., Mkn 501, lES 1959H-650). 
We now consider if this can also be the case for this flare of 
PKS 2155-304. 



ment of Thomson condition for X-ray photons {v^) form a closed 
system of three equations in the unknown variables 5 (beaming 
factor), B (magnetic field), and R (size) of the emitting region 
(see e.g.. iKatarzvnski et al.ll2005l Appendix A). These values 
can then be compared with the size of the emitting region ob- 
tained from variability timescales {R < ctvarS{l + z)~^), and 
with the condition that the cooling times of the electrons emit- 
ting at the peak {'jpeak) is equal to or shorter t han the variabil- 
ity tim escales (namely, jcooi < Ipeak', see e.g.. lTavecchio et akl 
fT998l) . 

With the parameters measured for this flare {vs ~ 50 eV, 
~ 0.5 TeV, is = 1-6 X 10*^ gj-g/^^ 1 x 10"*^ erg/s), 

the Thomson limit for photons in the middle of the Chandra 
passband {v^ — 1 keV) requires that b ~800, B ~0.6 mG 
and R ~lxl0^^ cm. This solution seems unrealistic: besides 
the quite extreme values of the beaming factor and severe effi- 
ciency problems, it is not consistent with the VHE spectra. The 
Thomson 7-ray spectrum should extend to 10 TeV with the same 
slope shown by the X-ray spectrum up to 1 keV, in contrast to 
observations. 

For electrons emitting at the SED peaks, instead, PKS 2155- 
304 is almost certainly in the Thomson regime. The Thomson 
regime is already satisfied for & k, 15, but the condition of trans- 
parency for 1-TeV ph otons (see e.g., iDondi & Ghiselliiiilll995l; 
iBege lman et an i2008h requires b > 30, adopting for reference 
the variability timescale of the main 7-ray flare it^ar ~1 hr). 
The Thomson regime for the SED peaks is also indicated by the 
absence of a significant change of (< 1.5x) in the presence 
of an inferred ~3 x shift of v^. Since in the KN regime the ratio 

Vcl Vs oc 7pg^j, oc Vs the change in Vc should have been ac- 
companied by an increase in Vs by 9x (vs. \/3x in Thomson), 
bringing it well within the Chandra band. 

Therefore, the SED peaks of PKS 2155-304 are in the 
Thomson regime, but a part of the Chandra passband is most 
likely affected by KN effects. With 5 > 100, consistency between 
the size of the emitting region estimated from both the observed 
Compton dominance and the variability timescale is obtained, 
yielding R «1 x 10^^ cm, B <5 mG and a Thomson condition 
at < 0.3 keV. This is also the result obtained by an analysis 
o n the B — 5 plan e follo wing the analytical approach described 
in lTavecchio et al.l ( Il998l) . 

However, this solution is not consistent with the cooling 
times of the peak electrons being equal to or shorter than the 
escape timescales. Because of the low values of the magnetic 
field and high beaming factors, in all the above cases the energy 
of the electrons cooling within the variability timescale {t' = 
t'varS/ (1 + z)), either by synchrotron or IC, is extremely high, 
namely 7coo; > 6 x 10^. This value is well above 7peafc ~9xl0^, 
meaning that the whole electron distribution up to the energies 
sampled by our observations did not have time to cool. 

In the one-zone frameset, therefore, the peaks of the SED 
cannot be explained by radiative cooling. This result is quite 
different from all previous estimates for the SSC parameters in 
PKS 2155-304, and is mainly due to the increase by a factor 
~10 in the separation between the synchrotron and Compton 
peak frequencies. 



8.1. One-zone SSC analysis 

8.1 .1 . Thomson condition and beaming factor 

In a single-zone, homogeneous SSC model, the knowledge of the 
frequencies of the synchrotron and Compton SED peaks and 
Vc, respectively), their luminosities (Lg and Lc), and the require- 



8.1 .2. Expanding blob 

An intriguing possibility is to explain the observed decrease 
of both synchrotron and IC emission as adiabatic cooling due 
to a rapid expansion of the emitting region from an initially 
very compact size ("explosion"), with the total number of par- 
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tides roughly constant. This hypothesis has two immedi- 
ate advantages: 1) it accounts for the initial high value of 
the Compton dominance as well as its rapid decrease (since 
in the Thomson limit Lc/Ls = U[.^^/Ub oc NR-'^{'y'^), 
where (7^) is the average over the electron distribution, see e.g. 
iGhisellini & Tavecchiol2008bh : 2) it can explain in a natural way 
the change of curvature of the VHE spectra through a change of 
internal transparency to 7-7 interactions. 

The initial region cannot be too compact, or otherwise it 
would not be transparent to 7-ray photons. However, in the case 
of PKS 2155-304, the steep and curved VHE spectra can allow 
for a moderate amount of possible internal absorption to be con- 
sidered (r-y-y ^^1-2, as obtained for example with S ^30-40), 
without requiring anomalous spectral shapes. The optical depth 
T increases with 7-ray photon energy as r(i/-y) = TQ{iy^ / vq)" , 
where a is the energy spectral index of the synchrotron target 
photons {F{i^) (X a = T — 1), and has a radial depen- 

dence cx R~^, for a constant A^. Therefore, internal absorption 
can explain the stronger curvature in the brightest state and a 
rapid expansion would decrease the optical depth making the 7- 
ray spectrum less curved at lower fluxes, as observed. In fact, the 
observed change of c urvature is opposite t o that expected from 
radiative cooling (e.g.. lMassaro et al ]l2008l) . though it could also 
be due to a change of the maximum electron energy. 

To obtain the cubic correlation, however, one must assume 
correlated variations of a second parameter, in order to shift 
the radiative output from the IC to the synchrotron channel. 
This can be achieved by v arying the magnetic field B (see e.g. 
ICoppi & Aharonianlll999l) : for a given particle injection, an in- 
crease of the magnetic field can lead to a reduction (cx B^^) 
of the flux of the inverse Compton 7-rays while increasing the 
synchrotron emission. 

Using the same formalism as in iKatarzviiski et al.l (|2005|) . 
it is possible to find the relation that yields Fj cx F^. The 
evolution of the radius R and magnetic field B can be param- 
eterized as R = Roit/toY and B = Bo(V^o)"" (see also 
lAtovan & Aharonianlll999l) . The adiabatic losses sustained by 
the particles during the expansion are taken into account. With 
the values for the spectral indices before and after the syn- 
chrotron peak of ai =0.7 (assumed for the optical spectrum) 
and a2 =1.4 (as measured in the X-ray band), respectively, a 
cubic correlation is obtained with m/r ~ —0.4 ( ±0.1 consid- 
ering ai between 0.5 to 0.9). Namely, the magnetic field must 
increase as the blob expands at a rate B cx i^^™/'" — 
and on the same timescales as the flux variations. This implies 
that the total energy of the magnetic field in the region must in- 
crease substantially, as Wb oc R'^^'^. Thus it should be created 
either locally (presumably by turbulent dy namo effects) or sup- 
plied from outside (see e.g., discussion in lAtovan & AharonianI 
[l999h . 

However, such strong and rapid amplification of the mag- 
netic field seems not consistent with the optical data. Rapid vari- 
ations of B within the emitting region are also bound to affect the 
synchrotron emission before the peak, produced by particles not 
yet cooled. With the previous parameters, the observed decay of 
a factor 2 in the X-ray band would cause a decrease of t he opti- 
cal flu x by 15% (compare F'^ (t) vs. F} (t) in Katarzvns ki et aT] 
l2005h . and on the same timescales as the X-ray and VHE varia- 
tions. This disagrees with the optical light curve, which is still 
slightly rising (by ~5%) and then remains basically constant 
during the overall X-ray A'HE decay (see Fig.|2|. While it might 
be possible to further trim this scenario to cancel out any vari- 
ation in the optical band, this would require an extreme fine- 
tuning. 
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Fig. 16. Possible SED of the new flaring component in 
PKS 2155-304, which would vary "up and down" beneath the 
persistent emission of PKS 2155-304 (see text). Blue triangles: 
SED of the flaring component at the flux maximum inside 
the X-ray/VHE simultaneous window. Red squares: expected 
level of the X-ray flux assuming a linear relation between X- 
ray and VHE flux variations. For reference, the same Chandra 
T300-Xmax spectrum and optical data shown in Fig. [TS] are 
plotted as well (highest black points). 



For all of these reasons, it is very unlikely that both of the 
peaks in the PKS 2155-304 SED during this night are produced 
by the synchrotron and SSC emission of one single zone. 

8.2. Two-SED scenario 

A simpler and possibly more reaUstic explanation is provided by 
the superposition of two SEDs, produced by two different emit- 
ting zones. The first is responsible for the usual "persistent" SED 
of PKS 2155-304, which, presumably, peaks in the UV and at 
a few tens GeV, but with rather low VHE fluxes. It is typically 
variable on longer timescales (^0.2-1 days), as shown for exam- 
ple by the long ASCA, BeppoSAX , and XMM-Newton observa- 
tions (Tanihata et al. 2001; Zhang et al. 2002, 2006). The second 
zone is responsible for the dramatic flaring activity of this night, 
and most likely of all the active period of July-August 2006. It 
needs to be much more compact and have higher bulk motion to 
account for the much faster variability timescales (0.05-1 hrs). 
The true variations of the X-ray emission from this second zone 
can therefore be as large as the 7-ray ones, but they are simply 
seen diluted in the "persistent" component, which have compa- 
rable or higher synchrotron fluxes. These variations are instead 
fully visible in the VHE band because there the contribution of 
the standard SED is at much lower fluxes. 

A two-zone scenario is rather common for explaining ma- 
jor flares in blazars. The main novelty of this event with respect 
to all previous HBL flares is that the bulk of the luminosity of 
the new component is now emitted in the Compton channel in- 
stead of the synchrotron channel. The synchrotron emission of 
this flare is not bright enough to "break through" the persistent 
emission and thus dominate both peaks of the SED. 

Within this scenario, one can try to estimate the SED of this 
new component at its brightest state during the simultaneous ob- 
servations, assuming that the flux of the persistent component is 
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roughly constant. This is shown in Fig. [16] In the 7-ray domain, 
the flux of the persistent component is supposed to be negligible 
with respect to the observed flare. The spectrum of the flaring re- 
gion is therefore given, to a first order, by the highest VHE state 
directly observed (the same as in Fig.fTSll. In the optical band, at 
the maximum of the flare the detected variation is ~15%, so the 
flux in the flaring zone should be 0. 15 of the observed flux. In the 
X-ray band, we have used the lowest-state spectrum observed at 
the end of the night as background file for the highest-state spec- 
trum. As a result, a spectrum with flux ^0.5-5 = 1.98x10^^° 
erg cm~^ s~^ is obtained, which is well fitted by a log-parabolic 
model with Tikcv = 2.17 ±0.06 and curvature 6 = 0.22 ±0.11. 
This gives the estimate for the synchrotron peak position around 
0.4 ke V, which is also consistent with the optical flux being much 
lower than the X-ray flux (see Fig.fTSll. 

This underlying SED yields comfortable values for the 
synchrotron-Compton modeling, and becomes consistent with 
the peak being due to radiative cooling. Such scenario also 
agrees with the spectral variations being larger in the VHE band 
than in the X-ray band (for the same diluting effect), and would 
naturally account for the behaviour of the optical data as well, 
if the optical rise is indeed associated with the VHE flare. The 
optical band would be reacting to the electron injection in the ris- 
ing phase but not in the decaying phase, due to the much longer 
cooling times. The variations at VHE can then be allowed to 
scale Hnearly or quadratically with the X-ray variations, as re- 
quired by the specific modeling. 

8.3. Diagnostic of tlie flaring region 

Although the ampUtude of the synchrotron variations of the new 
component is unknown, important constraints can still be de- 
rived on its properties in both cases of a quadratic or linear rela- 
tionship between X-ray and 7-ray fluxes. 

With the new parameters, a quadratic decay obtained through 
radiative cooling is readily possible with a pure SSC scenario 
because the peak of the synchrotron emission falls now directly 
within the Chandra passband. Indeed a solution in the Thomson 
regime can be obtained with the comfortable values of S > 30, 
B ~1 Gauss, and •jpeak — 3 x 10^. However, the high Compton 
dominance requires the emitting region to be very compact, with 
i? ~ 3 — 5 X 10^^ cm. The flare therefore would originate in 
a region whose size is of the same order as the Schwarzschild 
radius of the putative central black hole, which is esti mated to 
be ^10^ Mq (see discussion in lAharonian et al.1l2007al) . 

If instead the X-ray flux scaled linearly with the VHE flux, 
it implies that the new component has a high Compton domi- 
nance (w 20) constant in time, during all flux variations. This 
behaviour points towards an origin of the 7-ray peak from exter- 
nal Compton emission rather than a pure SSC mechanism. 

This is indeed expected in scenarios with a strong ra- 
diative interplay between different pa rts of the jet, such as 
the "N eedle/Jet" model proposed by iGhisellini & Tavecchiol 
(l2008al) : a compact and fast region moving throughout a larger 
jet sees the dense target field produced by the latter (which is 
thought to be responsible for the persistent SED), boosting sig- 
nificantly its Compton emission. Unlike the way it is envisaged 
by the authors, however, the SED of such needle should also 
contribute significantly in the X-ray band, in order to explain 
the small-amplitude but correlated variability. Alternatively, the 
dominant seed photons could come from regions further down 
inside the jet, where the fl ow has already strongly deceler- 
ated (i.e., at the VLBI s cale, iGeorganopoulos & Kazanasll2003l : 
iPiner & Edwardsl |2004|) . or by circumnuclear radiation fields 



if the flaring region is close to the central engine. In the lat- 
ter case, 7-7 interactions with these photon fields could leave 
a distinct abs orption feature imprinte d onto the broad-band 7- 
ray spectrum dAharonian et al.ll200 8a). which the combined ob- 
servations with Fermi-L AT and lACT might be able to reveal 
dCostamante et al.ll2007l) . More complex models with several in- 
teracting zones or jet stratification have also been developed, and 
rec ently appUed to the variability of PKS 2155-30 4 in 2006 (see 
e.g. iBoutelier et alj|2008l: iKatarzvriski et al.ll2008[) . They seems 
capa ble of reproducing som e features of the present observations 
(e.g.. lBoutelier et al.ll2008h . 

8.4. New mode of flaring in IHBL 

On a more general basis, it is interesting to compare this 
flare of PKS 2155-304 with the other major flaring events ob- 
served i n both X-ray and V HE bands, namel y from Mkn42I 
in 20 00 (iFossati e t alJl20Qg|), Mk n501 in 1997 ( iTavecchio et al.l 
llOOlt iKra wczvn ski et alj I2OO2I) and lES 1959H-650 in 2002 
dKrawczvnski et al]|2004l) . 

There are many common traits: they all display a pre-flare, 
"persistent" SED with the synchrotron emission peaking in the 
UV/soft-X-ray band, a steep X-ray spectrum that hardens during 
the flare pivoting around the UV/soft-X-ray band, and a flare 
luminosity approximately one order of magnitude higher than 
the typical source luminosity. The mechanisms and properties of 
the flare injection thus seem common. It is the radiative output 
that now differs significantly. 

In the previous events, such a high luminosity was emitted 
mostly through the synchrotron process, leading to a dramatic 
shift of the peak position in the SED according to the new peak 
frequency of the emerging component. The typical Compton lu- 
minosity, even at the flare maximum, has always been equal to 
or less than the synchrotron power (using the same EEL model 
for all sources). In this event, instead, the bulk of the flare lu- 
minosity is emitted through the Compton channel, yielding only 
minor modifications of the overall synchrotron emission. 

A bimodality therefore seems to emerge in the mode of flar- 
ing for HBL: either synchrotron dominated or Compton domi- 
nated, with the most extreme example possibly being provided 
by the "orphan flare" event of lES 1959H-650, which likewise 
occurred during one single night (June 4, 2002). 

It is intriguing to note that this difference might simply de- 
pend on the environment, namely on the location of the flaring 
zone with respect to the region responsible for the persistent 
SED: if the new injection/flare is taking place far away, there is 
little radiative interplay between the two zones, leading to a typ- 
ical SSC-type flare. When the flare occurs close to it, or close to 
the black hole where external fields are more intense, the out- 
come is external Compton-dominated flares. More campaigns 
targeted on the intra-night variability are needed to address this 
issue, but the diagnostic potential on the jet structure and the 
location of the "7-ray zone" is very promising. 



9. Summary and conclusions 

A full night of simultaneous, uninterrupted observations in the 
optical. X-ray and VHE bands was performed during an excep- 
tionally bright state of PKS 2155-304 in July 2006. A sampling 
of both light curves and spectra with unprecedented detail is ob- 
tained. For the first time among HBL, a high Compton domi- 
nance is observed, with a peak luminosity reaching ^10**^ erg/s. 
The variations in the X-ray and VHE bands are confirmed to be 
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highly correlated, both in flux and spectrally, but they follow a 
cubic relationship during the decay phase. Homogeneous one- 
zone SSC scenarios do not provide a consistent explanation for 
both the observed SED and the variability behaviour, indicating 
that a single particle population cannot be responsible for both 
peaks of the SED during this night. 

We have interpreted the data as the emergence of a new 
component in the SED, strongly Compton-dominated and thus 
without enough luminosity in the synchrotron channel to over- 
come the "persistent" emission of PKS 2155-304. This new 
component must be either very compact - of the order of the 
Schwarzschild radius of the putative central black hole (~10^ 
Mq)- if the emission mechanism is pure SSC, or external- 
Compton dominated, as would be expected in models with a 
strong radiative interplay between different parts of the jet. 

The richness and quality of the data obtained from this ex- 
ceptional campaign provide a fundamental testbed for a time- 
dependent treatment of the emission scenarios, which is needed 
to fully extract th e information on the physi c al conditions in the 
source (see e.g., Coppi & Aharonianlll999l: iKrawczvnski et all 
[2002; Boutelier et al. 2008). This is beyond the scope of this pa- 
per, and therefore left to future studies, but it is a necessary step 
forward with respect to more common parametric fits, which 
simply postulate the prompt electron spectrum and treat different 
states as uncorrected emissions. 

From an observational point of view, the results from this 
night raise the bar for the requirements of future multiwave- 
length campaigns on HBL, at least dur ing major flaring events . 
With flares as rapid as a few minutes jAharonian et alj|2007al) 
and flux changes of 1 order of magnitude in less than an hour, 
sufficiently long, continuous and strictly simultaneous observa- 
tions are mandatory for a meaningful analysis. Delays, gaps, 
or mismatches in the simultaneous coverage, even as short as 
30 minutes, can be fatal for the diagnostic potential. The mul- 
tiwavelength coverage should also care to sample the SED on 
both sides of the emission peaks: as shown in the present case, 
without the optical light curve some scenarios could not be ex- 
cluded. In this respect, Fermi-LAiT will provide the missing piece 
of information on the 7-ray hump for many objects, though typi- 
cally on longer timescales. The all-sky monitor provided by this 
observatory will be of great importance to catch the brightest 
Compton-dominated flares among HBL, triggering X-ray and 
lACT observations and thus allowing the study of this appar- 
ent flare bimodality (synchrotron or Compton-dominated) on a 
larger sample. 

The flaring activity of PKS 2155-304 in July 2006 demon- 
strates in a compelling way that the fast, intra-night variability 
is of fundamental importance to gain insights into the accelera- 
tion mechanism, emission processes, and location and environ- 
ment of the emitting regions, which are simply lost or smeared 
out with a sparser, long-term sampling. The final goal is to 
achieve full multiwavelength coverage of large-amplitude and 
ultra-rapid events such as those observed on the night of July 
27-28. 
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